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ABSTRACT 

We present a comprehensive analysis of holography for the bubbling solutions of Lin- 
Lunin-Maldacena. These solutions are uniquely determined by a coloring of a 2-plane, which 
was argued to correspond to the phase space of free fermions. We show that in general this 
phase space distribution does not determine fully the 1/2 BPS state of = 4 SYM that the 
gravitational solution is dual to, but it does determine it enough so that vevs of all single 
trace 1 /2 BPS operators in that state are uniquely determined to leading order in the large 
A^ limit. These are precisely the vevs encoded in the asymptotics of the LLM solutions. We 
extract these vevs for operators up to dimension 4 using holographic renormalization and 
KK holography and show exact agreement with the field theory expressions. 
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1 Introduction 



Supersymmetric supergravity solutions play an important role in developing our under- 
standing of holographic dualities. Such solutions provide valuable examples where one can 
carry out detailed computations and, using non-renormalization properties, make quanti- 
tative tests of gravity /gauge dualities at and away from conformal fixed points. Given the 
precise holographic dictionary available in such cases, one may also understand in detail 
how the spacetime is reconstructed from gauge theory data. One would hope to take from 
these examples generally applicable methods and principles, along with insight into the 
inner workings of holography. 

Supergravity solutions that asymptote to AdS^ x describe either a deformation of AT = 
4 SYM or the theory in a non-trivial state. The most supersymmetric non-trivial vacua of 
N = A SYM theory preserve 16 supersymmetries. In this context it is interesting to consider 
the SYM theory on both Minkowski spacetime R^^'^^ , in which case we have N = 4= on the 
Coulomb branch, and the theory on RxS^. These two cases are equivalent in the conformal 
vacuum since the two backgrounds are mapped to each other by a Weyl transformation but 
differ on a generic half supersymmetric state (which spontaneously breaks the conformal 
invariance of the = 4 SYM theory). Of course, the two theories are still related in the 
decompactification limit of S^. 

Using standard D-branc physics, one expects that the holographic dual of = 4 SYM 
on the Coulomb branch is the near-horizon limit of multi-center D3 brane solutions [1]. 
The Coulomb branch may be parametrized by the vevs of chiral primary operators and the 
gravity/gauge theory duality together with non-renormalization theorems imply that the 
vevs computed at weak coupling are non-renormalized and must therefore also be repro- 
duced by the holographic computation. In [2], building on [3, 4], we indeed succeeded in 
extracting these vevs from a generic multi-center solution, showing exact agreement with 
field theory. This provides a highly non-trivial test of the correspondence away from the 
conformal point - an infinite number of vevs was quantitatively matched - and also illus- 
trates the maturity of holographic methods as it shows that one can go beyond qualitative 
matching, performing precise quantitative computations. 

Supergravity solutions corresponding to AT = 4 SYM on RxS^ were recently constructed 
in [5]. The solutions of [5] preserve an i? x 50(4) x iS'0(4) bosonic symmetry and 16 
supersymmetries. These "bubbling solutions" are uniquely determined by a coloring of the 
2-plane into black and white regions. Based on earlier work relating 1/2 BPS states to free 
fermions [6, 7], this distribution was argued to map to the phase space distribution of free 
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fermions and supporting evidence was provided in [8, 9, 10, 11]; see also [12, 13, 14, 15, 16, 
17]. 

It is often stated in the literature that the gauge theory dual of the bubbling solution is 
the matrix model associated with free fermions. One of the aims of this work is to understand 
to which extent this assertion is valid by applying standard AdS/CFT methods. That is, we 
will address the question of whether the matrix model captures the entire vacuum structure. 
In the AdS/CFT correspondence the asymptotics of the supergravity solution encode QFT 
data. In particular the vacuum structure of the dual QFT can be extracted from the near 
boundary asymptotics of the solution. In the first part of this paper we will use holographic 
renormalization [18] and KK holography [4] in order to extract the vevs from the solutions 
of [5]. Any proposal for the field theory dual must reproduce these results^. 

Let us now discuss the QFT side. By the operator-state correspondence, all 1/2 BPS 
states of A'' = 4 SYM can be obtained by acting with 1/2 BPS operators on the conformal 
vacuum. The operators constructed from the 6 scalars X"^ of A = 4 SYM lie in the 
(0, /,0) representation of the S'C/(4) R-symmetry (see, for example, the review [19]). Up to 
a L'"(3) C SU (4) rotation, every such operator can be represented holomorphically using a 
single complex combination of the scalars Z = + iX'^ . Thus gauge invariant operators 
built from these scalars preserve an 5'0(4) part of the SU{4) R-symmetry and have a 
definite SO{2) charge j under rotations in the X^—X^ plane. A convenient basis for these 
operators is the Schur polynomial basis [6] and an arbitrary half BPS state |$) preserving 
SO {A) R symmetry can be written as a superposition of states 

m = ^aRXR{z)\n) = o^\n) (1.1) 

R 

for suitable complex coefficients or, where xr(^) is the Schur polynomial associated with 
the^ U{N) representation R and is the conformal vacuum. The representation R may 
be labeled by a Young tableau and the associated Schur polynomial Xr{^) has degree equal 
to the number of boxes n and in general involves both single and multi-trace contributions. 
Thus the operator 0$ is equal to a sum of terms each of which has dimension equal to 
charge, A = j = n, for any n> 1. It follows that in order to specify the theory we need to 
supply the coefficients ur and any gravitational dual should encode these coefficients. 

When the field theory is formulated on Rx one may reduce over the to obtain a 

^Here we assume that the 16 supercharges protect the vevs from acquiring quantum corrections, as in 

the case of = 4 SYM on R'^^-^K 

^We are really interested in SU{N) gauge theory but the difference between L''(A'') and SU{N) is sub- 
leading in the large A'^ limit. 
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one-dimensional model involving an infinite number of fields (KK modes). Given the large 
amount of supersymmetry, however, one might anticipate that the vacuum structure, i.e. 
the coefHcients aR in (1.1), may be encoded, at least to leading order in the large N limit, 
in the truncation of the reduction to only the s-mode of the complex scalar Z. We will 
take this as a working assumption in this paper. It would be interesting to investigate the 
validity of this assertion in general. We should emphasize that one should be very cautious 
about using properties of this matrix model which are subleading in (to infer properties of 
the dual spacetimes etc.) since these are likely to be different from the true 1/N corrections 
of = 4SYM onRxS'^. 

Standard arguments map the matrix model to free fcrmions, whose phase space dis- 
tribution is meant to correspond to the coloring of the 2-planc that determines the LLM 
solution. In the fermion picture, we have N free fermions with the ground state being the 
completely filled Dirac sea. This state corresponds to AdS^ x . Excited states are in one 
to one correspondence with the Schur polynomials, the excitation numbers being directly 
determined from the the length of the rows of the associated Young tableaux. A generic 
excited state in the free fermion picture is then in direct correspondence with the 1/2 BPS 
state (1.1) of AT = 4 SYM. The crucial question is then: does a phase space distribution for 
the fermions uniquely determine the state |$)? We show that this is not the case, i.e. the 
phase space distribution does not determine all coefficients a^, but it does determine the 
state enough so that the vevs of all (single trace) chiral primaries in this state are uniquely 
determined! These are precisely the vevs encoded in the asymptotics of the supergravity 
solution. 

The results of the holographic computation show that the LLM solutions (generically) 
encode vevs of all SO{4) singlet operators; such operators can be labeled by their SO{2) 
charge j. So to check the correspondence one should compute these vevs in the field theory. 
Here we face the first obstacle. While maximally charged operators (those whose SO{2) 
charge is equal in magnitude to the dimension) involve only the Z field and so can be 
implemented in the matrix model [10], all other operators involve all six scalars and thus 
appear to involve fields not included in the matrix model. 

To see how one can deal with this issue, recall that the 1-point function of an operator 
O in the state |<I>) is equivalent to the 3-point function in the conformal vacuum between 
O, the operator that creates |$) from \Q) and its conjugate, 

(o)$ = {n\oloo^\n). (1.2) 

Suppose this correlator is computed in free field theory. Since 0$ is constructed only from 
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Z the 3-point function receives contributions only from part of O that contains Z. Thus for 
the free field computation of the 1-point functions in (1.2) one may set to zero all fields but 
Z (and Z) in the chiral primary operators. We emphasize, however, that this truncation 
would in general give incorrect answers if used for different computations, e.g. two point 
functions. If the free field computation were to be renormalized, then fields apart from Z 
would of course contribute in loops. However, three-point functions of single trace chiral 
primary operators of A?" = 4 SYM are known not to renormalize [20] and it is believed that 
three point functions of protected multi-trace operators are similarly non-renormalized [21]. 
We indeed find that the vevs computed using free field results for multi trace operators 
do agree with those extracted holographically, thus confirming the expectation of non- 
renormalization. 

The truncated operators can therefore be implemented in the matrix model. We do this 
explicitly for all operators up to dimension four (whose vevs we also extract from the gravity 
solutions). In particular, we show that each of these operators can be expressed as linear 
combinations of bilinear s of fermion creation and annihilation operators. The coefficients 
in the linear combinations are fixed by demanding that the operators have zero expectation 
values at the conformal vacuum, 3-point functions with single trace operators are correctly 
reproduced, and the vev of the operators have the correct limit in the decompactification 
limit of the . In this limit the phase space distribution maps to the distribution of 
eigenvalues of the scalar s in the Coulomb branch oi N = A, SYM. Having implemented the 
operators in the matrix model it is then straightforward to compute their vevs in a general 
state 1$) and we find exact agreement with the holographic computations! 

This paper is organized as follows. In the next section we summarize how to extract 
holographic data from asymptotically AdS^ x solutions. The resulting expressions for 
the holographic vevs in terms of supergravity field asymptotics are applicable not just to the 
bubbling solutions of interest in this paper, but to more general 1/4 and 1/8 BPS bubbling 
solutions. In section 3 we review the LLM solutions and extract the vevs of all maximally 
charged operators and of all operators with any charge up to dimension four. In section 4 
we discuss the dual description of the bubbling solutions. We show what information about 
the state is captured by the distribution and hence the gravity solution; we reproduce the 
holographic vevs in section 5 and we explicitly match certain specific symmetric distributions 
with 1/2 BPS states in section 6. In section 7 we discuss our results. 

Appendices A and B review relevant properties of spherical harmonics and of scalar 
chiral primary operators in = 4 SYM whilst appendix C discusses the large N scaling of 
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three point functions. Appendix D is rather tangential to the focus of the paper: we discuss 
the KiUing spinors for the LLM supergravity solutions. These were discussed in [5] but 
only half of them were correctly identified and they are missing local phase factors which 
(drop out of the fermion bilinears used to construct the supergravity solution but which) 
are needed to solve the Killing spinor equations. 

2 Extracting holographic data 

In this section we will give a self-contained summary of the method of Kaluza-Klein holog- 
raphy, developed in [4], which allows the computation of all 1-point functions from any 
asymptotically AdSp x Xg supergravity solution. 

The basic steps in this method are the following. First one expresses the deviation of 

the supergravity solution from AdSp x Xq in terms of the eomplete basis of harmonics of the 
compact manifold Xq] let the expansion coefficients be denoted collectively as ip-^ . Now one 
forms gauge invariant combinations of these fluctuations, il^'^ , that satisfy field equations 
which can be expanded perturbatively in the number of fluctuation fields. Schematically 
these field equations may be written 

Cl4^ = CIJK^P^i^'^ + CIJKc4>^i>^^p'^ + ■■■, (2.1) 

where >Cxi - x„ is an appropriate differential operator. Since Cx^-i„ involves derivatives, the 
set of field equations cannot generically be integrated into an action. However, one can 
always define p-dimensional fields '^-^ by a non-linear Kaluza-Klein reduction map of the 
fields tp-^: 

= + IC^jK,^^^^ + • • • , (2.2) 

where /Cj-^ contains appropriate derivatives. The reduction map is such that the fields 
do satisfy field equations which can be integrated into an action. Given this p-dimensional 
action, it is then straightforward to obtain the one point functions of operators in terms 
of the asymptotics of the fields , using the well-developed techniques of holographic 
renormalization [22, 23, 24, 25, 26, 27, 28, 29]; for a review, see [18]. We will now give the 
details of each step in the case of interest. 

Let us consider any asymptotically AdS^ x solution of type IIB, which involves only 
the metric and 5-form field strength. (It is straightforward to include all other fields of type 
IIB, but unnecessary for our application here to the LLM bubbling solutions.) The IIB 
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SUGRA field equations^ for the metric and 5-form field strength are given by: 



Rmn = qFmpqrsFn^^^^ , 
These equations admit an AdS^ x 5*^ solution 



dsl = — + ^dx?^+de^ + sm^0dnl + cos^ed(f)'^ 



(2.3) 



(2.4) 



TpO 



'^IJ.vpaTi 



Fabcde — ^abcde- 



where (^, z^) and {a,b) denote AdS^ and indices respectively; M,N,... are lOd indices 
whilst X denotes AdS coordinates and y denotes coordinates. We will consider here 
solutions that are deformations of AdS^ x such that 



gMN = gMN + ^MN, 

Fmnpqr = Fmnpqr + Jmnpqr- 
These fluctuations can be expanded in spherical harmonics as: 



(2.5) 



and 



h(ab){x,y) 
K{x,y) 



fnupari^j y) 
fa^ivpcrix, y) 
fabfii/pi-^i y) 
fabcfj,u{Xj y) 
fabcdni-^j y) 
fabcde{x,y) 



Y.i4>ll^{x)Y^:^^{y) + ct/^{x)D^i^{y) + </>g)(x)D(„A)y^Hy)); 

Y,^'^{x)Y'^{y\ (2.6) 



= Y.^h%p-^x)DaY'^{y) + AD^^h['^^^{x)Yt{y)); 

= E(3^[Mfti] (^)^M (y) - 26;;i,(x)D[,y,f (y)); 

= Y.^2D^A^x)eabc''^DdYt{y) + 3b'^^ {x)D[,Y^'^\y)y, 

= Y.{D^bll^{x)eabc/D,Y'^{y) + (A^« - 4)b'; {x)eabcd'Yj^ (y)) 

= Y.h'^'^p)K''eabcdeY''{y); 



(2.7) 



Numerical constants in these expressions arc inserted so as to match with the conventions 
of [31]. Parentheses denote a symmetric tracclcss combination (i.e. ^(af,) = l/2{Aab + 
^ba) ~ l/Sfi'ab^a)- ' ^/^i ^(ab) ^[ab]" denote Scalar, vector and tensor harmonics whilst 
A^i and A^s are the eigenvalues of the scalar and vector harmonics under (minus) the 
d'Alembertian. The subscripts t, v and s denote whether the field is associated with tensor, 



^The field strength differs by a factor of 4 from the conventions in [30]. 
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vector or scalar harmonics respectively, whilst the superscript of the harmonic label 
derives from the number of components n of the harmonic. Relevant properties of the 
spherical harmonics are summarized in appendix A. 

In what follows it will be useful to label perturbations by both the degree k of the 
associated harmonic and by the degeneracy of such harmonics. For example, tt'^^ will 
denote the fluctuations associated with degree k scalar harmonics with I labeling the SO (6) 
quantum numbers. 

2.1 Gauge invariant quantities 

When computing the spectrum it is useful to impose the de Donder-Lorentz gauge choice, 
as in [31], which imposes the following conditions on the metric fluctuations 



along with analogous conditions on the five-form fluctuations. These gauge conditions 
remove terms involving gradients of spherical harmonics. 

As discussed in [4], it is often the case that the natural choice of coordinates for the 
asymptotic expansion takes the fluctuations outside the de Bonder gauge. Indeed, we will 
find here that there is a distinguished coordinate choice which is outside de Bonder gauge. 
This issue may be dealt with using gauge invariant combinations of the fluctuations; these 
were derived up to quadratic order in the fluctuations in [4] . For the purposes of this paper 
we will need only certain combinations which are gauge invariant at linear order, namely: 



Note also that /i^^ is a deformation of the background metric and it indeed transforms as 
a metric. 

2.2 The spectrum 

In this subsection we review the relevant parts of the spectrum of fluctuations about AdS^ x 
computed in [31]. As discussed in detail in [4], one can relax the de Bonder gauge fixing 
condition used in [31] by replacing all fields by the gauge invariant (hatted) versions given 
in the previous section. 



= D'^hai, = 0, 



(2.8) 




(2.9) 
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The scalars relevant here satisfy the fohowing Unearized equations 

Os''^' = k{k-4:)s''^\ k>2, 

□£^=^1 = (A; + 4)(A; + 8)Fi, k>0, (2.10) 
where we introduce the combinations 

with inverse relations b''^^ = -l^^^i + it^^^ = lOks''^^ + 10{k + 4)**^^!. The fields are 
dual to scalar chiral primary operators. 
The relevant vector combinations are 

= iBi;f,-^ik + 3)b'J^y, (2.12) 

with the corresponding masses being 

m2(a*) = (A:2 - 1); m\c'') = {k + 3){k + 5). (2.13) 

Thus the k = 1 modes of are massless and are dual to the R symmetry currents. 
The combination of lOd fields that satisfies the 5d linearized Einstein equation is 

= (hlu + Ig^'y, (2.14) 

the shift by tt^ follows from the Weyl transformation required to bring the 5d action into 
the Einstein frame. 



2.3 Kaluza-Klein reduction 

We are now ready to give the non- linear Kaluza-Klein map. This map relates the coefficients 
appearing in asymptotic expansion of the ten dimensional solution to the coefficients in the 
asymptotic expansion of the five dimensional solution obtained from the following action: 

/V-2 /• / 1 -I ^ 1 \ 

S = ^ y d^'xVG \--{R - 2A) + -G^^'^Y.^d^S^'d^S'^' + Y^F^^.F^n (2-15) 

where A is the cosmological constant and F^j^i, = d^Ajf^ — d^Ajf. In this expression we 
have kept only the fields that are dual to operators up to dimension 4, i.e. the metric G^i, 
which is dual to the stress energy tensor, the massless vector field Ajf which is dual to 
the R current and the scalars S'^^ dual which are dual to chiral primaries of dimension A;. 



10 



The reduction worked out in [4] included further low lying KK modes. However only the 
terms listed in (2.15) contribute to the formulae for the 1-point functions of the operators 
up to dimension 4. In particular, the potential terms for the scalars (derived in [4]) do not 
contribute. 

The relation between the fields appearing in this action and the 10 dimensional fields 
reads, 

S''^ = w{s'')s''^, k = 2,3; w{s'') = ^/8k{k - l){k + 2)z{k)/{k + 1), (2.16) 

S^i = ^fs^I-^^^(^83s^■r§^K^JJ^^g2JJJ^.^2K^. (2.17) 
5 \ 27z(4) J 

G^. = 9% + h% + L^.; (2.18) 

4" = (2-19) 

where z{k) is given in (A. 2) and a4j^2J,2K is the triple overlap between scalar harmonics 
(A. 3). In comparing these and subsequent formulae with those given in [4] and [2] note that 
the formulae in these papers were for the case of 5*0(2) x S'0(4) singlet harmonics for which 
(^40, 20,20/^(4) = 3/(2-v/5)- Here the hatted fields arc the gauge invariant combinations. For 
computing the vevs of operators up to dimension 4, one needs in general gauge invariant 
combinations quadratic in the fields and these are given in [4]. 

For the bubbling solutions discussed in this paper, however, linear gauge invariant com- 
binations suffice because, as explained later, the solution is in the De Bonder gauge to 
leading order and deviates from it afterwards. In particular, to the order we work, one can 
remove all hats from all formulae in (2.16)-(2.19), except for s'^^ for which one must use the 
linear gauge invariant combination in (2.11)-(2.9). 

The Kaluza-Klein relations in (2.16)-(2.19) result in an asymptotic expansion of the 5d 
fields to sufficiently high order so that the vevs of dual operators can be extracted, i.e. we 
obtain all terms up to order z'^ for the fields^ S^'', all terms up to order z^ for A^"- and all 
terms up to order for G^i^, where z is the Fefferman-Graham radial coordinate (see next 
section) . 

The reduction of gauge fields was not discussed in [4] but can be determined from the 

results of [32] for the quadratic action. Again, non-linear corrections to this reduction 

formula will not be needed in what follows since they will not affect the vevs of the R 

"^In [2] the KK map for S'^' was computed up to terms of order z'^ but these higher order terms are not 
needed here. 
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symmetry currents. The normalization of the 5d gauge fields is such that the corresponding 
R symmetry currents have the standard normalization, that is, their two point functions 
are given by [33] 

{Jt{x,)J^{x2)) = ^^^"'m, - d,dj)j^^—-^, (2.20) 
where 4d coordinates are labelled by and (a, b) label the 50(6) indices. 

2.4 Holographic 1-point functions 

The final step is to use the method of holographic renormalization to extract the vevs from 
the asymptotics of the 5d fields. This is by now a standard procedure except that here one 
needs to include additional terms to accommodate extremal couplings (see section 5.4 of 
[4]). The relation between field asymptotics and vevs is most transparent in Hamiltonian 
variables where the radius plays the role of time. The 1-point functions are then related 
to the radial canonical momenta of the bulk fields, which are expressed as (non-linear) 
functions of the field asymptotics [28, 29]. 

So to obtain the vevs one should use the following steps: 

1) expand the deviation of the ten dimensional solution from AdS^ x in harmonics of 
and in the radial direction; 

2) obtain the asymtpotic expansion of the five dimensional fields, using (2.16)-(2.19)); 

3) use the expressions for the exact 1-point functions in terms of the asymptotics of the 5d 
fields to obtain the vevs. 

This is the route followed in [4]. Since we know the Kaluza-Klein map in general one can 
process the 1-point functions to express them directly in terms of the asymptotic coefficients 

of the ten dimensional fields. Then the final formulae can be used directly after step 1) 
without having to know or use the Kaluza-Klein map. Such formulae were presented in [2] 
for a subset of operators and we present the formulae for all operators up to dimension 4 
here. 

Let us first summarize the expressions for the holographic 1-point functions in terms of 
the asymptotics of the 5d fields. The near-boundary expansion of the bulk metric G^u, the 
gauge field and scalar fields where k is the dimension of the dual operator, take the 
form 

dz'^ I 

= ^ + ^ {G(o)ij{x) + z^G(2)ij{x) + 2^(G(4)ij(a;) + logz'^H(^^)ij{x))) dx^dx^; 

Ai{x,z) = A(o)i(x) -Fz^A(2)i(x) H ; 

S\x,z) = z''(logz^Sl^ix) + Sl^{x) + ■■■); 
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S\x,z) = z(^-'=)sfo)(x) + --- + ^%,_4)(x) + ---, k>2. (2.21) 

In these expressions the boundary fields G(o)ij, ^(o)i, S^o)' "^fo) pai'ametrize the Dirichlet 
boundary conditions and are also the field theory sources for the QFT stress energy ten- 
sor, the R symmetry current and operators of dimension 2 and k, respectively. The gauge 
field is in radial axial gauge, = 0. The near-boundary analysis determines all coeffi- 
cients in these expansions except the ones corresponding to the normalizable modes, namely 

Consider first the scalar operators O52/ and Ogsi, where / labels their degeneracy. We 
will be interested in SO{A) singlet operators which can be labeled by their 50(2) charge m, 
but we will express the holographic relations in a more generally applicable way. For these 
operators the holographic relations are [4]: 

/V2 

{Os^') = 2^(4'))' (2.22) 

where 7r|^^ indicates the part of the canonical momentum of the field S''^ that scales with 
weight k. The relevant part of the canonical momenta can be expressed in terms of the 
asymptotic expansion of the 5d fields as follows^ 



TT, 



(2fc-4) 



(2fc-4)[5^^]fc (2.23) 



where the notation [A]if indicates the coefficient of the term in A and z is the Fefferman- 
Graham radial coordinate. The relation (2.23) holds for k ^ 2; when k = 2 one should 
replaces the factor {2k — 4) by 2. 

As discussed in some detail in [4] the vevs of the scalar operators of dimension four also 
involve quadratic terms; these are necessary to accommodate extremal couplings. Thus the 
vevs in this case are 

= £ (4') + 2^^414')0 (2-24) 

Using the Kaluza-Klein map in (2.16-(2.17) we now express these vevs directly in terms of 
the coefficients that appear in the lOd solution. 



^^'-'^ = £^[^''1^' (^5-) = £^[^'']3; (2-25; 



®Note that we chose not to include in the definition of tt the prefactor of Ar^/27r^, in constrast to the 
conventions of [4, 2]. 
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where z{k) is the normahzation of the degree k spherical harmonics, defined in (A. 2) and 
o-Ai,2j,2K is the triple overlap between scalar harmonics (A. 3). The expression for {Ogu) 
can be further simplified for solutions in which has vev (rather than source) behavior, 
such as those under consideration in this paper. In such cases, the asymptotics of (2.21) 
imply that 

[{D^s'''D^'s'% = [z''d,s^-'d,s^% = ^[s^'s^^ (2.26) 

and thus we obtain 

Next consider the stress energy tensor; its vev can be obtained by analyzing the coupled 

system of the metric and the scalar fields S^^ . (The other 5d fields fall off too fast to 

contribute to the stress energy tensor.) The part of the 5d action involving the metric and 

one field is same as the sector of gauged supergravity analyzed in [26, 27], where was 

called The result for the stress energy tensor can thus be carried over from these works, 

with S"^^ . Thus one gets 

jY2 / 1 ^ _ 1 

{Tij) = -^{G(^i)ij + -{Sl^'^Sl}^-^)G{Q)ij + -^^ (2.28) 

11 3 \ 

-2(^^(2))^ + 4<^(2)ijTrG(2) + 2-f^(4)ij ) > 

where the summation over / is implicit. 

We now want to rewrite this expression in terms of the coefficients of the ten-dimensional 
fields. To do this we should use (2.18) to express the 5d coefficients in terms of the coefficients 
of the ten-dimensional fields. To present a formula that is universally applicable wc should 
start from a universal form of the lOd metric. One such form is the Fefferman-Graham 
form, which can always be reached by a coordinate transformation. Note however that the 
expression (2.28) presumes that the bd metric G^i, is in the Fefferman-Graham form, but 
because Lzz,Lzi,n^ are generically non-zero, the Wd and 5d metrics cannot simultaneously 
be in the Fefferman-Graham form. This is not a problem since after the reduction to 
5c? we can always find the appropriate coordinate transformation that brings Gfj,i, to the 
Fefferman-Graham form. 

By definition the lOd metric is {g° + hP) but it is more useful to consider instead the 
following combination, 

9nv = 9% + (2.29) 
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and use as a starting point this metric in the Feffer man- Graham form, 

dz^ 1 

(2.30) 

The reason is that the z dependence of 7r° it is not a priori known, so had we started from 
g° + we would not be able to work out explicitly the subsequent transformation that 
brings G^i, to the Fefferman-Graham form. (Of course, in any given example one can always 
find this transformation, but here we aim to give formulae that apply universally.) 
Using 

J-'zz = S(o)S(o) ; J^zi= ^*(0)'^j*(0) "I ! 

and (2.18) we find in particular that 

= ^(1 + — z^s^Q-jS^o)), Gzi = --z^s'^^Q-^disf^^ (2.32) 
so the metric is indeed not in Fefferman-Graham form. The coordinate transformation 

z = z'il - y's%f%^ + •••), = + ^"c.~s%)di~s%) + • • • , (2.33) 

with a a STiitablc numerical constant, brings the metric to the Feffrcman- Graham form with 
the following coefficients: 

G{G) = 9(0), G(2) = 5(2), G{4) = 5(4) - ^^ioffo)9{o), ^(4) = ^(4)- (2.34) 
Inserting these expressions in (2.28) we get 



/ 2 

i'^ij) = ^ ( 5(4)ij - 9(^(0) 5(0) )5(o)ij (2-35) 
+^[Tr5(2) - 0^^9{2)f]9(o)ij - ^i9{2))ij + \9(2)ijT^g(2) + ^^(4)ijj , 



where we emphasise that g(k)ij the coefficients in the Fefferman-Graham expansion of 
the lOd metric in (2.29). . 

Let us finally consider the R symmetry currents; from the results of [27] their vevs are 



where again we rewrote the R symmetry current in terms of ten-dimensional fields to give 
the second equality in (2.36). 
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Before leaving this section, let us comment on the wider applicability of the highlighted 
expressions for the holographic vevs, (2.25), (2.35) and (2.36). In this paper we will analyse 
in detail the LLM bubbling solutions, which preserve an i? x SO{4) x SO{4) symmetry 
group, and are associated with 1/2 BPS states of A/" = 4 SYM on i? x S^. 

However 1/4 and 1/8 BPS states on RxS^ which are built from operators involving only 
the six scalars of A/" = 4 also induce vevs only for the R-currents, the stress energy tensor 
and the scalar chiral primaries. Thus the expressions for the holographic vevs given here 
can be used to extract such data from putative dual geometries, of the type constructed 
in [34]. It would be straightforward to derive corresponding expressions for more general 
asymptotically AdS^ x solutions, which involve more supergravity fields, such as the 
Janus solutions recently derived in [35]. Note in particular that the expression given here 
for the stress energy tensor (2.35) provides a rigorous way to extract the mass (including 
the Casimir term) from the ten-dimensional solution. 



3 Bubbling solutions 

The LLM bubbling solutions are^ 



ydyVi 

F5 

F 
F 
Bt 
dB 



-h-^{dt + Vidx'f + h^idy"^ + dxidx') + ye^d^l + ye'^dQl; 
2ycosh(G); z 



i tanh G; 



(3.1) 



€ijdjz\ y{diVj - djVi) 



F^t^dx^" A dx" A dOa + F^^^^dx^" A dx"" A d^s; 
dBt A {dt + V)+ BtdV + dB- 

dBt A {dt + V) + BtdV + dB; 



z + 

y2 



dB = -\y^ Hd 



where i = 1,2 and *3 is the Hodge dual on the parameterized by {y,xi,X2). The 
solutions are characterized by a harmonic function on six dimensions, with sources on an 
R^. That is, 

(3.2) 



z{x]_,X2,y) 



y' 



Jr2 ((x - X')2 +?/^)^ ' 



Note that we use the notation z with two completely diflterent meanings; as the function z defined in 
(3.1) and also as the Fefferman-Graham radial coordinate, (2.21). The meaning of z should be clear from 
the context. 
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where regularity of the solution requires that z{x'i,X2,0) takes the values The two- 
dimensional vector Vi can be written as 

e, , f z(x\,x'n,Q)(xj — x', ) dx\ dxn 

In polar coordinates on this can be written as 

r I' z{r' , (p' ,0){r — r' cos{(/) — (p'))r'dr'd(f)' 
~^ ^ ~^ Jr2 {{x-x'Y+y^Y ' ^ ' 



Vr 



z{r\<i)',Q)sin{(t)-(t)'){r'fdr'dct)' 

' ~ - 'r2 2/2)2 • 



3.1 AdS^ X solution 

The AdS5 x solution is obtained by taking sources for z on a disk of radius ro- Then 



^"=-h\ :y ... 1 (3-5) 



where r is a polar coordinate on i?^ such that xi = r cos (j) and X2 = r sin ^. Introducing 
the following coordinate change on the B? parameterized by y,xi,X2 



y = Rcos0 = RcosO; r = Rsm0 = \/ R'^ + r^sinO; (P = (j)-t, (3.6) 

gives 



2+ (^2+^2^0s2^)- ^^-n 

The coordinate shift (3.6) changes the flat metric on R^ 

dsl = dR^ + R^{dP + sin^ ed^"^) (3.8) 

to the following metric: 

2 +d9^) + (R^ + ^o) sin^ - dtf- (3-9) 

The other functions in the metric take the values 

(/i-2)o = ToHR^ + rlcos'^ey, {ye^y = rocos^ 9; (3.10) 

= - rp2'i'f (z/e-«r = -- 

(it^ + cos"' (9) ro 

The superscript A" denotes that these are the background AdS^ x functions, about which 
we will expand. Substituting these values into the metric gives 

/ jp2 \ 

= ro -{R^ + l)df + — ^ + R^dnl + {dO^ + sin^ Odcf)^ + cos^ 9dnl) , (3.11) 



(i?2 + 1) 
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where R = R/tq. This is indeed the metric on AdS^ x with global coordinates on AdS^ 
and curvature radius henceforth ro will be set to one. 

The five form field strength can be obtained in the following way. The two form dB is 
given by 

{dB)e^ = -{dB)0t = -\R^{R^ + l)cos^e sin OdR^; (3.12) 
{dB)R^ = -{dB)Rt = lR^cos^O sin ede^; 



where 



and for AdS^ x 



^ = y-\z + \), (3.13) 



^" R^{R^ + cos^ey ^^'"^^^ 
The two form dB is similarly given by 

{dB)e<p = -{dB)et = -^cosesin9{R^ + 1) - iR'^iR^ + l)cos^0sin9dR^; 
{dB)R^ = -{dB)Rt = -^Rsin^e+lR^ cos^ 9 sinOde^; (3.15) 

Substituting into the expression for the five form then gives the following expression for 
AdS5 X S^: 

F^R = R^; F^^ = cos^ 9 sin 9, (3.16) 

as expected. 

3.2 Asymptotic expansion 

Now let us consider more general solutions which arc asymptotic to AdS^ x 5^. The field 
theory data will be extracted from their asymptotic expansions around the AdS^ x 
boundary. This expansion can be economically expressed as follows. Let the solution be 
expressed in terms of the harmonic function $(xi, X2, y) with 

$ = $° + A$, (3.17) 

where $° is the harmonic function of the AdS^ x background about which we perturb. 
A$ can be expressed as 
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where Az(xi,X2,0) = {z{x'i,X2,0) — z°{x'i,X2,0)). Now note that $ (and hence A$) is a 
scalar harmonic function on which preserves SO (4) rotational symmetry. The asymp- 
totics can thus be expressed as 

A$(i?, e, ^) = ^(A$)fe„^^M, (3.19) 

k,m 

where Y^{6, (p) are normalized 50(4) singlet spherical harmonics of degree k with m label- 
ing their 50(2) charge; the properties of such harmonics are discussed in appendix A. By 
the addition theorem the coefficients in this expansion are given by [2] 

{A^)km = 2'=(fe + l)7r-i^^ Az(x;,4,0)(C-...,^x^i •••x*'0dxid4, (3.20) 

where C'^,,,^^ are 50(4) invariant symmetric traceless tensors on of rank k which are in 
one to one correspondence with the <S'0(4) singlet spherical harmonics. In particular 

(A$)2o = 4\/37r-i / Az{x'^,x'2,^){r'fdr'd(l)\ (3.21) 

where the explicit representation of the S'0(4) x 50(2) singlet tensor is used. In deriving 
this and subsequent expressions, one needs to integrate the explicit form of the normalized 
traceless tensor on R^, C^,,,^^, over the source distribution in i?^. Since only the coordinates 
{x^,x'^) are non-zero, terms in C^,,,-^, for which i 7^ 1,2 do not contribute. 

Note that the expansion (3.19) begins at A; = 2. There is no A; = term, since the 
leading asymptotics are those of AdS^ x and k = 1 terms are unphysical since they can 
always be removed by choosing the origin of the coordinate system to be at the centre of 
mass. The centre of mass conditions imply that 



z{r',4>', Q)r'e^''t>'{r'dr'd4>') = [ Az{r', 4>' , 0)r'e±''^'(r'dr'4') = 0. (3.22) 

JR2 



'R2 

The harmonic function (3.19) is expressed in the usual coordinates on R^, but to perturb 
relative to AdS^ x the function needs to expressed in terms of the coordinates [R^O^cj)). 
However, for {R,R) 3> 1, this change in coordinates changes the form of (3.19) only at 
k > 4, that is, 

A$(i?,^,^,t) = {A^h J'^'^^^t {A^)s J'"^^^;^t ~ + (3-23) 

{{A^)4mYr{0, <t>-t) + {A^)2mr'{0, ^ - t)) + ■ ■ ■ 

where the functions /"*(0, (t> — t) are given by 

f{e,(t^-t) = ftY^ + f2Y^{e) + ftY!{e)- (3.24) 
f^\e,^-t) = ft%^\e,cp-t) + ft^Y^\e,<p-t). 
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Note that Yo = 1. These coefficients are obtained by expanding A^{R,6,(f)) using 

1 1 / sin2 e \ ~ / cos2 e \ , , 

and then projecting back onto the basis of spherical harmonics. In what follows we will 
need only the following coefficients explicitly 

fo=0; /4° = -^; ft' = -'-^- (3-26) 

The functions appearing in the metric can be expressed in terms of A$ as follows 

ye^ = cos^ ^(1 + i?2(i^2 ^ pQg2 - cos^ ^(i?^ + cos^ 0)A$)-^ = cos^ ^(1 + a); 

ye-° = R^{1 + B?{B? + cos^ ^)A$)-5 (1 - cos^ e{R^ + cos^ 0)A$)5 = R^{l + (3)- 
h-\R^ + cos^ e)-^ = (1 + (i?^ - cos^ e)A^ - i?2 cos^ e{R^ + cos^ 0)2a$2)-2 = (1 + 7); 

+ cos^ ^) = (1 + (i?^ - cos^ e)A^ - R^ cos^ e{R^ + cos^ 0)2A$2)5 = {l + 5). 

Note that the leading order terms in (a, /? 7, (5) are of order ; to extract vevs of operators 
of dimension four and less it will be sufficient to expand these quantities up to order 
Then 

a = {\R^A^ + R^cos^eA^-\R^{A^f + ■■■)■ (3.27) 

8 

/? = (-^i?^A$-i?2cOs2M$+|i?8(A$)2 + ...); 

8 

7 = {-\R*A^ + lR\A^f + ■■■)■ 

o 

5 = aR^A^-\R^{A^f + ■■■). 
8 

Now consider the vector Vi: one needs to expand the perturbation of the vector written 
in polar coordinates (3.4) from the background value. Using the centre of mass conditions 
(3.22) one finds that (AF) = {V - V°) is given by 

AF^ = sin^ M$ + ^ (y''^\A^)22 + y''-2(A$)2(_2) + (A$)2o(2F2>o + _LyO^^ 

sin e (e2^^A$22 - e-2^^A$2(-2)) ■ (3-28) 



AVr = - ' 



Changing coordinates to {R, 9, (p) gives 



AV^ = -AVt = AV^; (3.29) 
AVr = sineAV;. + ---; AVe = RcosOAVr + ■ ■ ■ . 
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It is convenient to use a gauge transformation to remove the radial component of this vector. 
Shifting AV — AV + dA is equivalent to shifting the coordinate t ^ t — A whilst leaving 
^ = (f) — t unchanged. Choosing 

^ = -ji^iY'''iO,^){A^)22 - r2'-2(A$)2(_2)) (3.30) 

results in 

AVi = !:*(^ + ...; AF, = ^£<^ + .., (3.31) 

with 

v4e,<P,t) = -R^smHA^-'-d^{Y^'^{A^)22-Y^'-\A^)2^_2)) (3.32) 
+1 ((A$)2o(2y2'° + -^y°)); 

vg{e, 4>, t) = -'-de (r2.2(A$)22 - y^^-\a^)2^_2)) . 

Note that in these epxressions the spherical harmonics are functions of {9,(f),t). 

3.3 Expansion of metric and five form 

The asymptotic expansion of the metric is given by 



ds'' = -dt\{R^ + + sm^e{j-S) -2-^)+ R\l + P)dnl + {1 + S) 



dH 



2 



-2dtde^ + 2dtd(j){{j - 5) sin2 - ^) (3.33) 
+(1 + S)de'^ + cos^ ^(1 + a)dnl + 2 sin^ ^#(^) 

sin2^#2(l + <^ + ^(<^-7) + 2^) + ---, 

where the terms retained are sufficient to extract vevs of operators of dimension four and 
less. That is, to extract the vevs of the scalar operators with dimension less than or equal 
to four one needs (vr, (/)(5)) to order To extract the vev of the R symmetry current 

one will need -B(^)^ with /i / i? to order For the vev of the stress energy tensor one 

needs (tt'^, /ijjj^) up to order l/R^. 

In actually extracting these fields there is considerable simplification relative to the 
discussions of [4]. Consider first the perturbations tangent to the sphere and note that 

hahdx^dx^ = (1 + a){de'^ + cos^ Od^l + sir? Odcj)^) + 0{^). (3.34) 

R 

This implies that to order the metric perturbation along the sphere is in de Bonder 

gauge, with only the trace tt non-vanishing. Compared to [4] where the fields ^(5) (zero 
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in de Bonder gauge) were already excited at order there is simplification. Gauge 

invariant combinations of fluctuations are needed first at order (or correspondingly 

when computing vevs of operators of dimension four) and moreover only linearly gauge 
invariant quantities will be needed at this order. One would only need to use gauge invariant 
quantities at non-linear order for computing vevs of operators with dimension greater than 
four. Thus one can immediately extract 

^km ^^km^_^ ( A$),^e-^-* + • • • ; A: = 2, 3. (3.35) 
To obtain (tt^, vr^, 0^^^ we split the sphere perturbation into its trace and traceless parts: 
TT = (^5<5 + i?2(3-2sin2^)A$ + ^ + ---^ ; (3.36) 
V) = (-^i?'A$(3-2sin2 0)-|^ + ...); 
^(^^) = sin2 0Qi?2A$(7sin2^-3) + |^ + ---^ ; 
/i(^a^/3) = cos^ Og^^^p Qi?2A$(2 - 3 sin^ e) - |^ + . 

with he^ as given in (3.33), x° are coordinates on and g^a^ii is a unit radius metric on 
S^. Projecting onto the basis of spherical harmonics gives 

^An, ^ ^-.m. (^(A*)4„ - ^^A_„„^(A#,„A*2,)) + 1^^, + ' ' ' I (3.37) 



i40 



^ -(A$)20 + 



= (-^(A$2nA$,(_„)) + •••), 

where here and in subsequent formulae fflmpg = 0'A.m,2p.2q- The expressions for D(^f^Dy^Y^"^ 
given in appendix A arc used in extracting the values of ^(s). For example, in the case of 
the R charged harmonics the relevant perturbations are: 

V) = ^e±2^M(sin^0-sin2^)(A$)2±2; 

^00) = ^e±2i(<^-^)sin2e(-sin4e + sin2e)(A$)2±2; (3.38) 



V) = :^e±2^(<A-*)sin3ecose(A$)2±2, 
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which comparing with (A.22) and (A.23) gives 

A similar analysis follows for the neutral harmonics. Thus the gauge invariant combinations 



are 



^'"^ = e-^'"* - ^^^«-np(A$2nA$2p) + 200^76^-* + •••); (3.40) 

Now consider the vector fields. The metric (non-zero) fluctuations can be expressed as 

hta = {Y''\A^h2 - y^'-^(A$)2(_2)) - -^(A$)2ori. (3.41) 

The physical vector fields arise from the projection of the h^a terms onto vector harmonics 
to give 5^5^^. The non-zero projection of h^a onto scalar harmonics takes the metric outside 
de Bonder gauge, but the resulting vectors BjK do not contribute to any gauge invariant 
quantities computed here. Thus the only relevant vector term is 

Finally let us consider the metric perturbation. The perturbation h^^^ receives contributions 
only from the first line in (3.33). Thus the metric + h'j^^ with h^^ = h^^ + ^Tr^y"^ is 
given by 

ds^ = -di2(^i22 + l-^-^(A$)20 + ^(A$2nA$2(-n))) (3-43) 

+ (^^(^-3^^^^-^^^(-"))j 

where summation over n = (— 2, 0, 2) is implicit. 

Next consider the five form field strength. To compute the vevs we need only the 
modes {bj^, 6^^^ in the expansion. This means that we need only expand /^^ = {F^a) and 
fe4> = {Fect> - F^^) giving 



fe4> = Sine cos^ e (^\R^{R^ + l)dRi^^ + -^{l-2,s\n^e) + s\necose 



2E? 



cos 6* 



fuB = sin0cos3 0(-ii?3^^ + ... ) ; (3.44) 

sm 6 



fn^ = sm9cos^e{lR%A^ + ---)- 
fet = sm9cos^e(^lR 



fet = sm9cos^e(lR^dRiA^) + 2a-:^^dea 
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Prom the fot term one gets 



fet = sinecos3^(-^-^(A$)2o + •••), (3.45) 



from which one can extract 

Combining this with the vector (3.42) extracted from the metric one finds that 

at =-^(A<I>)2o; 4 = 0. (3.47) 

This is the anticipated result since the massive vector should not be excited at this order. 
Prom the fotf, terms one finds 

bf"^ = -in^^-'"'\^^)km; k = 2,3 (3.48) 



and therefore the gauge invariant quantities are 

Putting together (3.40) and (3.49) gives 

= ""^'"'^i^^^^)^- + (^^^'"VS - 19 j4^(4) »^"p(^^2nA$2p))44 + •••)• k = 2, 

= "l92^I(4)""™*'*-"^'(^^2'^^^2^)- 

Note that there are no terms in the gauge invariant fields t^"^. This is a computational 
check: using [4] these fields satisfy the field equations 

(□ - 96)t^"* = 32z(4)-ia^nps2"s2p^ (3.5I) 

and thus at order l/R^ can only receive contributions quadratic in A$2n- The terms 
are linear in A$2n and thus cannot contribute to the fields t^"^ at this order. 

3.4 Holographic vevs 

Given the asymptotic expansions of the relevant fields we can extract the values for the vcvs 
using the formulae from section 2.4. The relation for the R symmetry current vev (2.36) 
along with (3.47) implies that 
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To apply the formula (2.35) for the vev of the stress energy tensor one must first bring the 
metric (3.43) into Fefferman-Graham form, by the coordinate change 



Then 

where ga/s is the metric on the unit radius 5"^. Using the explicit form for (A$)2o from 
(3.21) and reinstating factors of a, the inverse radius of the S^, gives 



(Jt) = 






(3.56) 


(Tn) = 




[^^ pir^ - W)rdTd<i^ = (r„), + a{.h) 





where {Tu)c is the Casimir on Rx and the density function p(r, (p) satisfies 

f p{r, (f))rdrd(t) = 1; p° = -^e{a - r). (3.57) 

We define 9{x) = 1 for .x > and 9{x) = otherwise. A general distribution is such that 
p{r, (f)) takes the value I/tto? in a region of the plane with area vra^, and is zero everywhere 
else. The corresponding mass E and R-charge J are given by integrating these expressions 
over the 5^, resulting in 



J = AT^a / p{r'^ - W)rdrd(l)- (3.58) 
E = N'' + / , P^''^ - h<'^ydrdci^ =Ec + aJ. 



These quantities have the expected behavior, namely J = for AdS with the Casimir 
energy Ec taking the expected value; the energy and angular momentum tend to zero in 
the limit of a large S*'^ and the BPS bound {E — Ec) = aJ is saturated. 

For the scalar operators (2.27) along with (3.50) implies the following result for the vevs: 

(k - 2) Uk-l) 



22'=(A; + 1) ^ ^ ^ 



80i?Vg^ = -^(A$)2±2; 80i?V^,^) = --^(A^>)2o, (3.59) 
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with (A; — 2) — 1 for A; = 2 and the scale a, the inverse radius of the S^, reinstated. For 
operators with \m\ = k (and k ^ 1) the vevs are therefore 

(C>5fc±fc) = ^^(fe-2)Vfe^e-*"** / {r^p)e^'^'^rdrd(t). (3.60) 



Recall that there is no /c = 1 operator in the SU{N) theory; the integral vanishes in this 
case because of the centre of mass condition (3.22). For the other operators with dimension 
less than four, one gets 





^/2iV2 
\/37r2 j 


' {r'^Az)rdrdcf); 


(3.61) 




V2N^ / 
V37r2 \ 


J p{r^ — ^o?)rdrd(j)^ ; 






9 


/ {r^p)e^''l'rdrd(t), 













where in the second expression for the neutral operator the explicit form of p° is used. For 

the opcra(;ors with dinwrnsion four, again rc^instating thc^ inA'cu'sc radius of S"^ mxc finds 



(O540) = 


/ Az(3r^ 4a\ydrdct); 

VStt^ Jr2 


(3.62) 








(C>54±2) = 







The general structure of the vevs is thus 



i(fc-|m|) 

{Osk^} = iV^e"'™'^* Y p{r^-'^U^^)e''''^rdrdct), (3.63) 

with certain coefficients ai . These vevs can also be written in the form 

{Ogkm) = >ffee-^"'"*(A$)fe^ + akma^{Os(k-2)m) + hmo" {O s(k-,)m) + ■■■ (3.64) 

where (A^^, ctkm) appropriate constants. In the a ^ limit only the first term survives, 
and as will discuss below one recovers the Coulomb branch result. Since k > \m\ the 
vevs of maximally operators only receive contributions from the first term. Our explicit 
computations go up to dimension four, but if one assumes this structure persists in the vevs 
of higher dimension operators then the result (3.60) holds for maximally charged operators 
of all dimension. We will find that the field theory result does indeed reproduce (3.60) for all 
k, thus verifying this hypothesis. By contrast the vevs of non-maximally charged operators 
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do receive other contributions and thus one needs to calculate explicitly the appropriate 
coefficients {akm,Pkm, ■■■)■ 

These expressions make manifest the Umiting behavior as a ^ and the theory passes 
to that of the Coulomb branch of = 4 on i?^'^. The R-charge and the energy as given in 
(3.56) vanish in this limit, as expected for supersymmetric vacua of A/' = 4 on R^'^. However, 
the scalar chiral primary vevs remain non-trivial for appropriate density functions p(r, (/>). 
Each density function describing a regular bubbling geometry consists of N droplets dj, 
such that p{r, 0) takes the value I/tto^ on the droplet, and the area of each droplet is tto^. 
Suppose the boundary of the droplet is described by r = r j + di{^), with rj constant and 
some suitable function ((?!)). Then the density function describing the droplet is 

PdM 0) = -^ein + di{(P) - r), (3.65) 

such that 

Pdi {r, (j))rdrd(l) = —. (3.66) 

Coulomb branch solutions arc then obtained in the limit that stays finite as a — > 0: the 
density function for each droplet behaves as 

Pd^ir, </.) ^ ^5{x^ - xDSix"" - xf), (3.67) 

satisfying (3.66). Here (xj,xj) describe the location of the droplet in the 1-2 plane, and in 
this limit each of the droplets is associated with an eigenvalue of the matrices {X^, X'^). 
Clearly in the a ^ limit the disc density function describing the conformal vacuum 
becomes a delta function localized at the origin, p{r,(j)) 6{x^)S{x'^). 

Now taking the a — > limit in the vevs of the scalar chiral primaries one gets 



/ 

Jd 



2'=/2(fe - 2)^ 1^ dx^dx^p{x\x^){C^^^.,^x'^ ■ --x^ 



7r2 V ft Jr2 

in exact agreement with the Coulomb branch vevs given in [2], restricting to an <S'0(4) 
invariant distribution. 



4 Dual description 

In this section we will consider half BPS states in = 4 SYM, and their relation to free 
fermions. We discuss the correspondence between an arbitrary half BPS state and a two- 
dimensional density distribution, which in turn is to be identified with the defining density 
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function of the bubbling supergravity solution. In particular, we show that the state is 
not completely determined by this density distribution, but the density distribution does 
determine uniquely the vevs of all single trace chiral primary operators. These in turn are 
precisely the information that is encoded in the asymptotics of the LLM solutions. Thus 
one would anticipate that the LLM solutions receive higher order corrections, involving 
information beyond the density function, which capture the dual state uniquely. 

4.1 Half BPS states in Af ^ 4 SYM 

There is a onc-to-one correspondence between half BPS 50(4) symmetric representations 
of A/' = 4 SYM and symmetric polynomials in the eigenvalues of a complex matrix Z or 
Schur polynomials. Here Z is one combination of the six Hermitian scalars X'^ of = 4 
SYM, given hj Z = + iX^. 

There are several choices of basis for the gauge invariant multi-trace polynomials of Z: 

1. The trace basis of products of traces of Z is an obvious gauge invariant basis. For the 
group U{N) the multitraces can be labelled by p{n) conjugacy classes of the permu- 
tation group Sn where p(n) is the number of partitions of n. Labeling representatives 
of different conjugacy classes of Sn by aj, the basis of multi-trace operators is given 
by TtiaiZy. 

Tr(a/Z)=yzf •••Zf . (4.1) 

jl—jn 

For SU {N) Z is traceless and one must therefore restrict to elements of Sn without 
1-cycles; the distinction between U(N) and SU{N) is however not important in the 
N ^ oo limit relevant here. 

2. The Schur polynomial basis is, in the case of U{N), a sum over these trace oper- 
ators, weighted by the characters of a in the representation R oi Sn, namely 

XRiZ) = -,Yl XRicrmaZ). (4.2) 
n! 

The representations R can be labeled by Young diagrams with n boxes, which corre- 
spond to partitions of n and there are thus p{n) Schur polynomials of degree n. An 
advantage of this basis is that the two-point functions are diagonal. Again modifica- 
tions are needed for the case of SU{N), but these give 1/N effects which will not be 
relevant here. 

Note that another useful basis is the dual basis, dual to the trace basis, but this will not 
play a role here. 
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An arbitrary half BPS state |$) preserving SO{A) R symmetry can therefore be written 
as a superposition of states 

1^) =Y.^RXR{Zm) =Y,biTj:{aiZ)\n), (4.3) 
R I 

for suitable (complex) coefficients a^i and bj, with \Q) being the conformal vacuum. Denot- 
ing by O"^ the set of gauge invariant operators, the vevs of these operators in the state |$) 
are given by 

{o^U = J2^>R'(^\(xR(^))^^^XR'{zm; (4.4) 

R,R' 

= Y.b}bj{n\{Tr{aiZ)yO^Tr{ajZ)\n). 

i,J 

A state is an eigenstate of the dilatation operator and of the R-symmetry (in the 1-2 
directions) with eigenvalue n if and only if the superposition involves only Sn- That is, only 
operators involving n fields Z are included in the superposition. 

It is also important to note that in the N ^ oo limit there arc considerable simplifi- 
cations in three point functions appearing in (4.4). Let us consider first computations in 
the trace basis, where the operators are normalized as C^iTr{a^Z) with n the dimension. 
The normalization factors C„i are such that the basis is orthonormal in the large N limit, 
namely 

C,^C,j^{n\{Tr{a^Z)yTr{afZ)m = Sjjd^^ + 0{1/N), (4.5) 

and the large N scaling of C^j is l/iV". 

Now consider the three point functions (4.4) in which the operators O'^ are single trace 
operators built from the six scalar fields X"^. These are clearly the relevant operators to 
compare with the holographic results. As discussed in [36] three point functions which are 
extremal, so that the conjugate operator has a dimension which is the sum of the dimensions 
of the other operators, and those which are non-extremal are known to have different large 
A*" behavior. Since the state |$) is a sum of terms each of which is maximally charged, 
i.e. it has j = A, it follows that the 3-point functions arc never extremal when O"^ is not 
maximally charged. Moreover the large N behavior depends on whether the other operators 
in the correlator are single or multi-trace. We discuss in appendix C the large N behavior 
of such correlators, and summarize here the relevant results: 

Non-extremal correlators: Non-extremal three point functions for which O-^ are 
(orthonormal) single trace operators scale as 1/A^ or smaller in the large A^ limit. Note 
that this assumes that the dimensions of all operators in the correlator are small compared 
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to N. In the case that the operator O"^ is neutral under SO{2) x SO{A) R symmetry 
the correlators behave as 1/N only for the diagonal terms, namely when cr" = ctj*. When 
the operator O"^ has a non-maximal SO (2) charge m the correlators involving single trace 
operators still behave as 1/N; correlators involving multi-trace operators are generically 
subleading in N, but for special cases can also behave as 1/N. 

Extremal correlators: Extremal three point functions scale as one or smaller in the 
large N limit. Correlators involving only single trace operators behave as whilst 
multi-trace correlators in which Tr{a^Z) = 0'^Tr{ay' Z) (and therefore a'j is necessarily 
multi-trace) are of order one. 

One can also rephrase these results in terms of the Schur polynomial basis, as discussed 
in [6]. As we will show in section 6.1, the vev of the S0{2) x 5*0(4) neutral operators in the 
state built from a given Schur polynomial of dimension n is independent of the choice of 
Schur polynomial. To leading order in it behaves as n/N. This result has an immediate 
corollary: consider geometries dual to different superpositions of the Schur polynomials, all 
of the same dimension n. Then symmetry implies only the neutral operators acquire vevs 
but these vevs differ only hy 1/N effects, so these geometries are not reliably distinguishable 
within supergravity. 

Now consider a superposition of states of different dimension (and thus R charge). In 
such a case SO{2) charged single trace operators acquire vevs, and the computation of 
the vevs of maximally charged operators necessarily involves extremal correlators. The N 
scalings of these vevs depend crucially on the specific Schur polynomials, or equivalently 
multi-trace operators, appearing in the superposition. Superpositions involving single trace 
operators will lead to vevs which are suppressed by 1/N relative to multi-trace superposi- 
tions, and thus these are immediately distinguishable. An explicit example illustrating this 
effect will be discussed in section 6.3. 

4.2 Relation to free fermions 

Consider irreducible representations of the symmetry group Sn- These may be characterized 
by a sequence of non-negative integers {A} = (Ai, • • • ,Xn) with Ai>A2>--->AAr>0 
and Ai = n. The sequence defines a Young tableau with the number of boxes in the 

ith row being A^ and the total number of boxes being n; let x^xyi^) be the corresponding 
Schur polynomial. In this section we will review the relationship between Schur polynomials 
and free fermions. 
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We introduce a second quantized free fermion field 

oo 

nz,z*,t) = Y,Cie-'^'^''^%iz,z*) (4.6) 

1=0 

where (Q, Cj) satisfy the anti-commutation relation {Ci, Cm} = Sim- Note that throughout 
this section we will set the inverse radius a of the to one; this sets the mass scale in the 
matrix model, and thus of the fermions, to one. The functions ^i{z, z*) are the orthonormal 
wavefunctions of the lowest Landau level, and are given by 



M^,^*) = ]/^z'e-''\ (4.7) 
The fermion field ^{z, z* ,t) satisfies the constraint that the total number of fermions be A^, 

/oo 
dzdz*^''{z,z*,t)-^{z,z*,t) = ^cjCi = N. (4.8) 

1=0 

The ground state is denoted and is given by 

\n) = Cj,_,Cj,_,---Clcl\0), (4.9) 

where |0) is the Fock vacuum defined by QjO) = for all /. We will denote by |<I>) a generic 
state containing N fermions; each such state can also be expressed as a superposition of 
Schur polynomials. The Schur polynomial X[x} (■^) corresponds to the state 

CI^,^_CI\0). (4.10) 

Now consider the expectation value of the density function defined as 

U{z,z*,t) = '^\z,z*,t)'^iz,z*,t). (4.11) 

In the conformal vacuum 

N-l 

{U{z,z*,t))n = Y.^i(z,z*m^,zl; (4.12) 



1=0 
N-l 



= E ^{zz*ye-^''* = 27r-'e-^''*EM-2{2zz*), 



=0 



where by definition 



For iV » 1, 



E^_,{y) = f2i- (4.13) 



e-yEN-i{y)^e{N-y) + f{y), (4.14) 
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where 0{x) = 1 for a; > and 9{x) = for x < 0. The function f{y) describes the smearing 

of the step function; /(y) has support only within a region around y = N oi width of order 
one and is such that f{y)<0 for y < N, f{y) > fov y > N with 

/ dyf{y) = 0; / dy\f{y)\ = 0(1). (4.15) 
Jo Jo 

Thus to leading order as ^ oo 

{U{z,z*))n = -9{N-2\zf). (4.16) 

TT 

To compare with the supergravity results one therefore needs 

^H; {U)a = 2p, (4.17) 

with identified with the supergravity coordinate r. In a generic state |$) the density 
function is given by 

{Uiz,z*,t))^ = J2^iiz,^*)^miz,z*)e'^'-"''^\cjCmh; (4.18) 

l,m 



l,m 



where we define 
The supergravity density function is related to this via 



Ut = {CjCmh- (4.19) 



{U{t = 0))$ = 2p. (4.20) 
4.3 Extracting the state from a distribution 

In this section we consider how to derive the specific superposition of Schur polynomials 
corresponding to a given distribution. Recall that the expectation value of the density 
function given in (4.18) is 

(?7(z,z*,t))$ = 5;a>;=(z,z*)cI.^(z,^*)e^('-™)*(C7;C'^)$. (4.21) 

l,m 

Thus by integrating a given density function with respect to a suitable basis of orthonormal 
polynomials one can extract the coefHcients (C/Cm)*- Note however that these expansion 
functions 'I>^<I>m are not orthogonal, when integrated over the plane with unit measure. 
So in practice it is actually more convenient to work with the density function in phase 
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space, u{p,q,t), which is naturally expanded in a useful basis of orthonormal functions. 
The explicit relationship between the density functions U{z,z*,t) and u{p,q,t) was given 
in [10]: 

U{z,z*,t) = J ^^e-^*^+^^*-^^* J dpdge-^(^+*f)+^*(«-^^')u(p,g,t). (4.22) 
Then, following [10] one finds that 

J dzdz*i-l)''izyz>'Lf = (e-i^^* J dpdge-^(«+^f)+^*(«-^P)n)^ ^ (4.23) 

Now to make manifest the behavior in the large N limit one should rescale these coordinates 
as in (4.17) so that 

, , , [n 

Fl = YyFl; Q=\j^x- p=^—y. (4.24) 

Retaining only the leading order terms in (4.23) as N ^ oo gives 

j cPw(w*yw^U = j d\r^+'=e^('=-^>u, (4.25) 

where x = r cos (j), y = r sin (p. As we have seen from the holographic computations, and 
will discuss below, all one point functions are expressed in terms of these integrals. Thus 
at leading order in N one can identify = r and the difference between the distributions 
{U, u) is not visible. Whilst it is more natural for the droplet distribution in the bulk solution 
to be identified with the phase space distribution, rather than the z space distribution, this 
is not distinguishable at leading order in N. 

In some calculations, such as that of one point functions which we will discuss below, 
it is more convenient to use the z space distribution, and exploit the simple form of its 
expansion in exponentials. For the current purpose it is rather more convenient to work 
with the phase space distribution, since this is expanded in a natural basis of orthonormal 
functions, the Laguerre polynomials. That is, the phase space distribution is given by [10]: 



m<n 

+ E \/^(-l)"x^^"^~"^e-^>^e^(--")^L--"(x)(C'Jn4)$, 



m>n 



where x = Nr and "^(x) is the Laguerre polynomial defined by 



n + a \ 



Kix)=Y.(-ir{ (4.27) 



p=o \ n-p 
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for which the orthogonahty relation is 

/ dxe-^x''L-{x)L^{x) = 5mn^-—r^, (4.28) 

Jo IT- 

for integral a. In the conformal vacuum one gets 

N-l ^ 

'iTpn = Y.{-re-2^LM (4.29) 

m=0 

where L^(x) = L^{x)- Using the identity 

dxe-2^Lm{x) = 2(-l)"^ (4.30) 

one can show that this satisfies the normalization condition f cPxpQ = 1. Moreover for 
large the distribution asymptotes as before to a disc, vrp^^ 9{1 — r). 

Using the orthogonality relation for the Laguerre polynomials one can now extract the 
{CLCn}^ via: 

{ClC^+p)^ = (-irl^l^^ I dcPe^P^dxx'^'e-^^'LMP^; (4-31) 

{cUpCmU = i-irl^^^^^ J d<l>e-'P'^dxx''/'e-^/'Kix)p^, 
where {m,p) > 0. 

Thus from the distribution one can extract the complete set of (C'mC'„)$. Let us 
now discuss whether knowledge of these is in principle sufficient to determine the state |$)^. 
Consider first the case where |$) has definite dimension n, so that 

l^)=E«{A}KW)- (4-32) 

{A} 

Normalization of the state implies that Yl{x} = 1- such a state {CliCp}^ is 

non-zero only for m = p and 

{CiiCm)'S> = E |«{A}P'^{A}m) (4-33) 
{A} 

where 5^x}m = 1 ifi^ the corresponding state contains a fermion at level m. Therefore 
one cannot extract the phases of a^x} from this information: the density function is not 
sufficient to completely determine the state. There is one exception to this: when precisely 
N of the {ClCm) $ are non-zero and equal to one, the corresponding state is necessarily 
a single Schur polynomial. In this case the summation in (4.32) collapses to one term, 
'^Related discussions appeared in the recent paper [38] . 
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and the overall phase of the state plays no role. Note that this is precisely the case that 
was discussed in [37], but for a general state of definite dimension the distribution is not 
sufficient to determine the state. To determine the phases in the general case one would 
need to know in addition 

iUcUW^m,)^, Y.mi = Y,mi. (4.34) 

i j i j 

It will be made manifest in the next section that {CmC.p)<^ determines the expectation values 
of single trace operators in the state, see (5.32), whilst (4.34) is related to the expectation 
values of multi-trace (neutral) operators. 

Note that the discussion so far has made no restriction on N. Even if one can determine 
the {CmCm)^ exactly one can still not determine the state. Of course when one takes the 
iV — oo limit and sharpens the distribution such that it gives a regular supergravity solution 
the situation will be worse. One will not be able to determine the coefficients {CmCm)^ 
exactly, and thence one can only determine the leading behavior in N of the |a{;^}.p. 

Now consider a general state |$) which does not have a definite dimension, so that 

l^)= E «n,{A}l";{A})- (4.35) 

n,{A} 

The neutral vevs {CmCm)<^ still do not determine the phases of the a^^^x). However, some 
phase information is obtained via the vevs of charged operators. That is, 

(C^+pCm)* = ^ al+p,{x}^an,{\}{n + p;{\}p\Cl^+pCrn\n-,{\}), (4.36) 

n,{A} 

= XI "n+P,{A}p"«,{A}' 
n,{A} 

where the Schur polynomial \n + p;{\}p) is precisely Cl^_^pCm\n;{\}) . That is, the state 
\n + p;{\}p) differs from |n;{A}) by only one fermion. So only a subset of the phase 
information is obtained; this is sufficient to determine the state when the superposition 
contains Schur polynomials such that each of which differs by only one fermion from at 
least one other polynomial in the superposition. If however the state contains at least one 
Schur polynomial which differs by two fermions or more from all other Schur polynomials 
in the superposition, then the phase of the coefficients of these terms cannot be determined 
without vevs of multi-trace operators. Thus for a general distribution one would again need 
vevs of multi-trace operators to determine the coefficients in |^>) uniquely. 

The supergravity solutions are constructed entirely out of the density function p$ and 
therefore contains information only about the expectation values of single trace operators. 
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To determine the state |$) one needs the expectation values of ah other operators, which 
are not determined by p^. Therefore, one would expect that the higher order corrections to 
the LLM bubbling solution, apart from correcting the distribution would also involve 
additional information so that the corrected solution captures the entire vacuum structure. 
This is in line with the fact that the IIB supersymmetry rules are expected to receive 
non-trivial higher derivative corrections. 

5 Computation of vevs 

In the previous section we set up the correspondence between a general half BPS state 
and a density distribution. The information extracted holographically is the vevs of chiral 
primary operators, and in this section we will discuss how these vevs may be computed in 
an arbitrary half BPS state. We compute explicitly vevs of all operators up to dimension 

four, and the vevs of maximally charged operators of arbitrary dimension, and find exact 
agreement with the holographic results. These results arc a detailed confirmation of the 
correspondence between LLM bubbling solutions and 1/2 BPS states, and moreover provide 
evidence that the vevs are not renormalized, as one might anticipate given the sixteen 
preserved supercharges. 

5.1 Energy and R-charge 

In a generic state |$) the energy and R-charge J relative to that of the conformal vacuum 
is 

7V-1 

{E-E,) = J = Y, m{ClCmh - E = E MClCm)^ - h^iN - 1), (5.1) 

m m=0 m 

where Ec = SAT^/IG is the Casimir energy on RxS^. This Casimir is clearly not reproduced 
correctly by the matrix model, since there are contributions from all KK modes on the 
of all SYM fields. Now note that 

/ d2^|^P(J7(z,z*,t))$ = ^^^|±^((7t,C'„)$. (5.2) 

m 

Thus 

m{ClCmU = [ d^z{2\zf - \){tj{z, z*,t))^, (5.3) 

m 

and hence 

{E-Ec) = J = j (fz{2\z\^ -\{N + l)){U{z,z\t))^; (5.4) 
= j dMH2-l(l + l))p, 
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which agrees with the holographic result (3.56), after taking the N ^ oo limit. 



5.2 Vevs of maximally charged operators 

We now consider how the vevs of single trace operators O^''' of dimension A and 5*0(2) 

charge k may be computed. Here we will use O^'^ to denote the operator in field theory, 
whose vevs we compute within free field theory, whilst OgAk refers to the corresponding 
operator whose vevs at strong coupling were computed holographically. 

Let us begin by computing the vcv of the maximally charged single trace scalar operator 
of dimension k, O^'^, in a generic state. The operators O^'^ are implemented as follows: 



1=0 

The factor X^^k is chosen such that the two point function satisfies the following normaliza- 
tion condition: 

{{0^^M)\(t^)0^-'^\t2)) =Ml5k,k,- (5.6) 
where A/fe is defined in (B.4). A/fe is the appropriate normalization for the two point functions 
extracted holographically. The normalization condition implies that 

^^^(/ + A;)! 1 {{N + k)\ N\ 



The corresponding integral representation of the operators is therefore {k ^ 1) 

O^'^'it) = Mk\k,k2-^ J dzdz*z''U{z,z*,ty, (5.8) 

TV , k k f 

= -—=y/k^{k -2)22 N-^ dzdz*z^U{z,z\t). 

TT y K J 

The expectation value of this operator in a generic state |$) is then given by 

(0*^''=(t))$ = -^^Vk^ik-2)2^N~^ f dzdz*z^{U{z,z\t))^. (5.9) 

The integral may be rewritten using using (4.17,4.20) as 

(C»'='*^(i))$ = ^y/k^{k - 2)e*'=* / d^we'^'t'^^kp^ (5_;L0) 

vr^V A; J 

in exact agreement with the holographic result (3.60). 

Consider the k = 1 operator. Here the distinction between U{N) and SU {N) becomes 
important: the vanishing of the trace in the latter means that there is no dimension one 
operator, O^'^ . This constraint can be incorporated here by restricting to configurations in 
which 

(CJi'i) = = J d^wwp, (5.11) 
which is indeed the condition imposed on the holographic distribution. 
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5.3 Other scalar chiral primciries 

Now let us consider the remaining scalar chiral primary operators O"^^, which are not 
maximally charged, A > \k\. The action of such an operator on the conformal vacuum 
or any other 1/2 BPS state |$) built from Schur polynomials creates a state which cannot 
be described in terms of Schur polynomials. The reason is that the other scalar operators 
contains not only the scalar fields (Z, Z) but also the remaining four A/" = 4 SYM scalar 
fields. The latter are not contained in the Schur polynomials, and thence not in the free 
fermion description. 

Suppose however one wishes to compute one point functions of these scalar operators in 
a state |$). To do so, following (4.4), one needs to know three point functions between such 
an operator and two maximally charged operators. Consider the computation of such a 
three point function in free field theory; at tree level the computation actually only involves 
the fields {Z,Z). Take for example a three point function such as 

{Tr{Z)'^{x)O^P'\y)Tr{Z)'^{z)), (5.12) 

for which the single trace SO (2) x 5*0(4) singlet operator has the structure 

O^f'O = aiTr {{ZZf + •••)+ aaTr {{ZZy'^R'^ + ...)+..., (5.13) 

where the ellipses within the trace denote cyclic permutations, (ai,a2,---) are constants 
and = Yli=ii-^i)'^ denotes collectively the other scalars of = 4 SYM. Since the 
latter have no propagators with the fields (Z, Z) and cannot be self-contracted, only the 
first term contributes in the three point function (5.13). 

Therefore, one would anticipate being able to implement the scalar operators with free 
fermions such that one can compute such one point functions. One would not however 
expect to be able to compute two point functions, or general higher point functions of such 
operators, using the free fermion description. 

5.3.1 Neutral operators 

Suppose one implements the dimension two neutral operator as 

©2,0 ^ ^^^2,0 CLCmim - (3). (5.14) 

m 

Then the normalization factor A2,o and the constant (3 should be fixed such that the one 
point function of this operator vanishes in the conformal vacuum; the three point function of 
the operator with charged operators gives the correct J\f = 4 results and the vev reduces to 
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the Coulomb branch result as the radius of the is increased. Imposing the first constraint, 
(C)2.0)^ = 0, implies that 

iV-l 

/3 = AT-i ^ m = i(Ar - 1). (5.15) 

m=0 

Note that this value for /3 implies that the operator O^''^ annihilates the conformal vacuum, 
= 0, and therefore the two point function for this operator also vanishes. This makes 
manifest the point made above, that one can only obtain the vevs of neutral operators from 
the matrix model. The corresponding integral representation of the operator is 

C>2,o ^ 2M2X2,oJ dzdz*{\z\^ - \N)U{z,z\t) (5.16) 

where subleading terms as N ^ oo are dropped. The normalization A2,o is fixed by taking 
the flat space limit: only the leading order term is retained, and comparison with the result 
(3.68) gives 

A.„ = ^. (5.17) 

We have also checked explicitly that this result is consistent with = 4 three point func- 
tions, involving charged operators. Using (4. 17), (4. 20) one can rewrite the vev as 

= ^ / ^'^d^l' - h)P^ (5-18) 
in exact agreement with the holographic result (3.61). 

Now let us apply the same techniques to obtain expressions for the vev of the dimension 
four neutral operator. For the operator C>^'° one gets 

= ^2^ ClCmi'-^ + hm + h). (5.19) 

Here the overall normalization is again fixed, so that the operator in the integral represen- 
tation gives the correct expression in the Coulomb branch limit. Imposing the vanishing of 

the vev in the conformal vacuum implies that 

2N^ + N^{Abi - 3) + iV(862 - 46i + 1) = 0. (5.20) 

Calculating the three point function with single trace charged operators of dimension two 
such that §2 = J\f2^0'^''^ gives 

where the ellipses denote terms which are subleading as AT — oo. Then solving these 
equations to leading order in N gives 

bi = -N; b2 = IN^. (5.22) 
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These values can be shown to also be consistent with other three point functions involving 
different charged operators. In integral representation this implies that 

= J dzdz*{3\zf - 2N\z\^ + \N^)U{z,z*,t), (5.23) 

where again only leading terms as AT — > oo are retained, and using (4.17) this gives 

(0^'°) = ^^J^ / (fw{3\w\^ - A\wf + l)p, (5.24) 

which agrees with the holographic result (3.61). 

5.3.2 Charged operators 

Now let us treat the charged operators in a similar fashion. For the dimension three operator 

O^'' = ^J dzdz*z{\z\^ + C3)Uiz, z*,t), (5.25) 
which implies in the fermion representation 

^ ^'7!^''* 5] ((m + 2 + 2c3)(m + 1)^^) cl^.Cm. (5.26) 

Now we use the three point function (B.8) to fix the coefficient C3; to leading order in N 
this gives 

cz = -\N, (5.27) 

and thus 

(03'i) = il^e'* / d^ww{\w\'^ - l)p. (5.28) 



7r2 



However the constraint (5.11) implies that the second term in (5.28) vanishes, and thus that 
the holographic result (3.61) is reproduced. 
For the dimension four operator 

= j dzdz*z\\z\^ + c,)U{z,z\t), (5.29) 

which implies in the fermion representation 

O''^ = e'^*-^ ^ ((m + 3 + 2c,W{m + 2){m + l)) Cl_,,Cm. (5.30) 



10A^2 

III 

Again we use a three point function (B.8) to fix the coefficient C4; to leading order in N 
this gives C4 = — ^N, and thus the vev is 

4^3 A'2 



iO''') = '-^e^^'l d'ww\\w\^ - l)p, (5.31) 
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in agreement with (3.61). 

Thus, to summarize, we have implemented all single trace operators of dimension A and 
SO {2) charge k quadratically in fermions as 

^ oo 
m=0 

with P^'^[m) fixed so as to give the correct normalization and three point functions of 
J\f = 4: SYM. Up to dimension four, it was sufficient to use three point functions involving 
only single trace operators, but for higher dimension operators one might also have to use 
additional three point functions involving multi-trace operators. Rewriting the vevs of these 
operators as integrals over the distribution gives the general form for the holographic vevs 
(3.63) and explicit agreement for all operators up to dimension four and maximally charged 
operators of all dimension. 

6 Correspondence between supergravity solutions and states 

In this section we explore how much one can deduce about the dual state from a given regular 
supergravity solution, using the information about the vevs of chiral primaries. We have 
already argued that even the exact distribution function does not in general determine the 
state uniquely, and in this section we will see how a given sharpened distribution (which gives 
a regular supergravity solution) can correspond to a number of distinct exact distributions. 

We will also note that non-singular supergravity solutions which break the SO {2) rota- 
tional symmetry are necessarily dual to infinite superpositions of Schur polynomials. Super- 
positions of a small number of Schur polynomials typically give rise to distributions which 
cannot be approximated by step functions and thus do not correspond to regular geometries. 
Thus the natural field theory bases for half BPS states, which use R charge cigenstates, are 
not the natural bases for describing regular bubbling geometries. 

We illustrate this point by considering a disc distribution with a ripple deformation of 
frequency n. Using the chiral primary vevs we argue that such a distribution is given by 
a coherent superposition of single trace operators. This identification follows very natu- 
rally from earlier discussion of quantum Hall liquids: area preserving deformations of a 
disc droplet are naturally described by coherent superpositions of fermionic excitations, or 
equivalently in terms of a collective chiral boson description. 
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6.1 Radially symmetric distributions 

In this section we consider half BPS states associated with superpositions of Schur polyno- 
mials of the same dimension. Such states are eigenstates of the dilatation and R charge, so 
only S0{2) neutral operators can acquire expectation values. The corresponding distribu- 
tions therefore preserve the rotational symmetry in the plane. 

A given Schur polynomial x^x^{^) corresponds according to (4.10) to a distribution 

N 

{U{z,Z*,t))n,X = Y.K+N-pi^'^*^t)^X,+N-p{z,Z*,ty, (6.1) 
p=l 

^ 2 |,|.^ (2|.p)^^+^-^ 

Given such a distribution is radially symmetric, only the stress energy tensor along with 
neutral operators acquire expectation values. The expectation value of the former is clearly 
independent of A, and depends only on n. One can now show that the vevs of neutral oper- 
ators with dimensions 2k <^ N also do not distinguish between different Schur polynomials, 
for A'^ — > GO. First note that 

Jdz\z\ {AUiz,z,t)U, = 2 (A, + iV-p)! 

= E (^(^ + k-p)-^{N- p))Xp + ■■■ = 2-^'A^-ifcn + • • • (6.2) 

p=i ^ ' 

where AUn^x = {Un,x — Ufi) and ip{x) is the Digamma function and ellipses denote terms 



which are subleading as A ^ cxd. Thus the leading term depends only on "^—i A 



p=i '^p 



and not the specific Schur polynomial. The vevs of neutral operators can be expressed in 
terms of such integrals as 

{0'''\,x = ^Y.dij <fz\z\'^'^-^^N^{AU{z,z\t))^,x, (6.3) 

for certain coefficients d/. (In the previous sections we gave the di explicitly for A; = 1,2.) 
Using (6.2), one finds that 

(0^'='°)n,A = = Y.d{l-^^-'\k - I), (6.4) 

regardless of the choice of A. Note that this behavior concurs with the explicit result for the 
vev in the state created by the single trace operator Tr(Z"), given in (B.6), and the latter 
result determines that 

\/2k: 

Ck = , (6.5) 

2'=-iV2^TT 

42 



Figure 1: Sharpened distributions describing Schur polynomials of the same dimension 
n <^ N cannot be distinguished by the corresponding vevs at leading order in N. 



Thus the vcvs of neutral operators in an R-charge eigenstate are at leading order in A'" 
independent of the specific choice of Schur polynomial superposition creating that state. 

Expressed in the coordinates appropriate for comparing with supergravity, the density 
distribution takes the form 



The latter implies that the excess energy relative to the conformal vacuum is n, indepen- 
dently of A. This density function is such that < 7rpn,x < 1 everywhere; however, as 
for the density function describing the conformal vacuum, (npn^x) does not take the values 
{0, 1} everywhere, and therefore the corresponding supergravity solution constructed from 
Pn,x would be singular. Just as for the conformal vacuum, though, the density function can 
be written as a sum of thcta functions plus correction terms describing the smearing which 
are sub leading as N oo. Retaining only the former leads to a non-singular supergravity 
solution. 

For example, suppose one considers the Schur polynomial for which Ai = n + 1 with 
Ap = otherwise; this corresponds to the state Clf_^_^CN-i\^)- The density function (6.6) 
in this case consists of a smoothed disc of radius one, along with a second peak localized 
around \w\'^ = {1 + ^^). Now suppose one sharpens distribution so as to get a regular 
supergravity solution, consisting of a disc plus an annulus: 




N 



(6.6) 



and has the properties 




(6.7) 



P = - ■ 

TT 



1 



0<\w 



(6.8) 
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P = 0: (1-1)<H^<(1 + ^); H^>(l + ^). 

The smearing of the disc and the annulus to obtain the exact density function is described 
by correction terms which are subleading as A?^ — ^ oo. 

As a second example, consider the Schur polynomial for which Ai = ni and A2 = ^2 = 
n + 1 — ni, corresponding to C'jy^^^j_iC'jy_j_^^_2C'Ar-iC'jv-2|^)- Assuming that ni and n2 
differ by a finite amount, the density function (6.6) consists of a smooth disc, along with 
two localized peaks. (If n\ and n2 are comparable, these two peaks merge, to look like 
the one.) One can then obtain a corresponding supergravity solution by taking the density 
distribution to consist of a disk plus two annuli: 

(.,(-_!)),„.,(,,-). 

Both configurations we have described have the same energy (n+ 1), and by the arguments 
above also have the same one point functions at leading order in N. 

More generally, for an arbitrary Schur polynomial one can obtain a corresponding super- 
gravity solution by sharpening the distribution into a set of annuli, as illustrated in Figure 1. 

However, there is clearly not a unique map from such a set of annuli to a given Schur polyno- 
mial: Schur polynomials which are very similar to each other, and superpositions of similar 
Schur polynomials, give density distributions which can only be distinguished at subleading 
order in N . That is, the associated sharpened supergravity distributions, which consist only 
of annuli, can be the same. To give an example, consider a specific superposition of Schur 
polynomials of the same dimension 

n 
fe=0 

with Yllt=o l^fcl^ ~ When only one coefficient is non-zero this collapses to the case 
discussed above of a single Schur polynomial, for which the sharpened distribution is a 
disc plus annulus. However, many other superpositions, such as those for which one is 
much greater than the rest, will give precisely the same sharpened distribution. Moreover, 
a typical Schur polynomial for which 1 ^ n ^ in which many of the \p are non-zero 
and different can give rise to a distribution which does not approximate a disk plus annuli: 
generically there are no strong peaks at \w\ > 1, as illustrated in Figure 2. By contrast a 
Schur polynomial for which the Ap are equal does give rise to a strong peak, see Figure 3. 
For more discussions on these issues, see [11], [37] along with the recent paper [38]. 
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Figure 2: Distribution for a typical Schur 
polynomial; there is no distinct peak. 
The figure shows A'' = 100, n = 30, with 
a random distribution of Xp. 



Figure 3: Distribution for a Schur poly- 
nomial in which the Xp are equal. The fig- 
ure shows N = 120, n = 30 with Xp = 2 
for p < 15. 



Note furthermore that a giant graviton, which is also a Schur polynomial of given di- 
mension, necessarily corresponds to a radially symmetric distribution and not a disc with 
an excited droplet illustrated in Figure 5, as suggested in some earlier papers. Indeed the 
maximal giant graviton, which has dimension N and is a Schur polynomial for which all 
Xp = 1, corresponds to a distribution which approximates a disc with a hole in the middle, 
as illustrated in Figure 4. 

In the above discussion we have focused on Schur polynomials with dimension less than 
N. A Schur polynomial with dimension greater than or equal to N can never be represented 

by a single trace component, since there is of course no single trace operator with such a 
dimension. Indeed if one considers the two distributions (6.8) and (6.9) but now with 
n > N, one finds that the vevs for the dimension four neutral operators are respectively 

(O540) = ^(2n + 3n2); (O540) = ^(2n + 3(n2 - 2n2ni)), (6.11) 

where n = n/N, hi = ni/N and n2 = n2/N with h = hi + h2 > 1- The leading order terms 
as N ^ 00 for (ni,n2) 3> 1 have been retained. These vevs are indeed distinguishable 
even as iV — 00; note that the vev for the single annulus is greater than that of the other 
distribution. Of course, following the general discussion given earlier, this distinguishability 
will still not be sufficient to determine the dual state |$) uniquely: the single trace operator 
vevs cannot provide enough information. Moreover, again a typical Schur polynomial will 
not give rise to a distribution with strong peaks which is well approximated by annuli. 

Consider a disc plus annulus configuration with fixed energy n. For given n one might 
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Figure 4: A maximal giant graviton corresponds to a disc with a small hole at the centre. 

think that the radius of the annulus \w\'^ = 1 + (3 can be taken to be arbitrarily large 
provided that its width is decreased accordingly. However, one has to take into the flux 
quantization condition: as discussed in [5] one can integrate the five-form flux over a sphere 
surrounding the annulus, and the flux must be quantized. This quantization requires that 
the total area of the annulus (and indeed of any isolated droplet) must be a multiple of 
tt/N. Thus the width of the annulus defined as 6 = \w\'^g^^ — \w\'^^j^ is a multiple of 1/N. 
The total energy of the configuration of the disc plus annulus, relative to the conformal 
vacuum, is given by n = NP + 1 and therefore for given n one gets (3 = {n — 1)/N, as 
in (6.8). This is the maximal radius for a given energy: multiple annuli or thicker annuli 
necessarily have lower radii. 

This same bound is also visible directly from the density distribution. First note that 
the function e~^x°'/(^^- ^oi large a has support only around x = a. Thus the distribution 
(6.1) extends to a maximum radius 

which is maximal for the Schur polynomial in which Ai = n and Aj = otherwise. This 
implies an upper bound on the magnitudes of vevs of the scalar operators, for a given energy. 

6.2 Non symmetric distributions 

A distribution which is not radially symmetric corresponds to a superposition of R-eigenstates. 
The generic state |$) can be written in terms of Schur polynomials as 

1^) = E«MX{A}I^^) ^Y.^n,x\n,m- (6.13) 

n,A n,\ 

with the Schur polynomials being orthonormal, and normalization of the state implies 
X^nAl^",Ap = 1. Now consider an operator Caj of dimension A, R-charge j. The vev 
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Figure 5: A disc plus droplet corresponds to a superposition of Schur polynomials of all 
dimension. 

of such an operator in this state is 

{OAjh= <+j,X'(^n,x{{n + j),{>^'}\OA,j\n,{X}). (6.14) 

n,\,\' 

The three point functions appearing in this sum arc non-extrcmal, except when A = j, i.e. 
the operator is maximally charged. As discussed previously, the leading order contribution 
to non-extremal correlators as A'^ ^ oo is independent of the specific choices of Schur 
polynomial. Thus, the leading order contributions to the vevs of non-maximally charged 
operators can be computed using 

|$)«^an|n), (6.15) 

n 

where |n) is any state of R-charge n, and the coefficients a„ are such that 

J2K,xf = \anf- (6.16) 

{A} 

However, the vevs of the maximally charged operators do depend on the specific Schur 
polynomials appearing in the state |$), because the corresponding three point functions are 
extremal. Thus these vevs can be used to distinguish between different distributions with 
the same a„, even at leading order in N. For example, consider the states 

1$) =ao|J^)+a„,{A}|n,A), (6.17) 

where |aoP + IfJn.lAjP = 1- To leading order in N, the vevs of the energy, R-charge and 
neutral operators arc independent of the specific choice of {A}. However, the vevs of the 
maximally charged operators with dimension n clearly do distinguish them. Let us compare 
the state created by the single trace operator, namely \n,X) where the single 

trace operator is defined in (5.5), and that associated with the Schur polynomial for which 
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Ai = n, namely |n, A) —>■ Clf_^^_iCN-i\^)- Then the expectation values of the single trace 
operator of charge — n are respectively 



1$) = ao\n) + anO'^'^in) : (0"'"")$ = iV^a^a^e"^"*; (6.18) 
1$) = ao\n) + anC},^^_^CN-i\n) : {C»"'-")$ = ^a*a„e-^"*, 

which differ by a factor of ^/n. 

We should note here, however, that such superpositions which involve only a small 
number of Schur polynomials do not generically give rise to smooth supergravity solutions. 
In the example just given, taking (ao,a„) to be real, the corresponding distribution takes 
the form 

p{w, 4>) = p{\w\) + p{\w\) cosincf)), (6.19) 

with the functions {p{\w\), p{\w\)) dependent on the specific Schur polynomial. In the case 
that |n. A) — > C']y^^_^C'Ar_i|ri) the functions are 

M\M) = Sd-iV-'-lwl^j+lwl^C-ye-l-l" (Kl'(f^+l°»l VC'-l)! l'°l'") • 

AriV+n/2-1 

' ' v^(iV-l)!(Af + n- 1)! ^ ^ 

Regularity of the supergravity solution requires that TTp{w, cj)) takes the values {0, 1} every- 
where. This condition can however never satisfied by a function of the form (6.19), with 
p(|i(;|) non-zero. One can see this easily as follows. Suppose the function p{w,<p) satisfies 
this condition at = 0, so that 

irpiw, 0) = p{\w\) + p{\w\) = {0, 1}, (6.21) 

for p{\w\) > 0. Then 

Trp{w, Scf)) - Trp{w, 0) = -^Tm^p{\w\)S(f)'^ + ■■■ (6.22) 

Regularity would require that the right hand side takes only the values { — 1, 0, 1} for all \w\ 
and 6<p, but this is not possible given that 6(p^ is a continuous function of 5(j). Smoothing the 
distribution so that TTp{w,(p) does take the values {0, 1} everywhere introduces additional 
terms (with small coefficients) into the state |$), and such terms may not be distinguishable 
at leading order in N. 

Conversely, a droplet distribution which is not radially symmetric but which gives rise 
to a regular supergravity solution is always associated with a superposition of an infinite 
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Figure 6: A disc with a ripple of frequency n corresponds to a superposition of Schur 
polynomials whose dimensions are multiples of n. The actual figure has n = 8,a = 0.1. 



number of Schur polynomials. Such a distribution can be written in terms of step functions 
whose arguments are the boundaries of the droplets: 

e{\wf-fm; f{cl>) = fo + m, (6.23) 

where /o is a zero mode whilst f{(f>) has no zero mode and can be expanded in Fourier 
modes. Now even when the droplet boundary f((p) contains only one frequency n, the 
distribution will contain all multiples of this frequency. One can see this by mode expanding 
the distribution, using (4.31), or by computing the multipole moments of the distribution 
as 

J #d|u;|2|u;|2'e^"^^p = (1 + Z)"^ J #e^"^^(/o + /(</'))^+^ (6.24) 

The latter integral is clearly non-zero when the frequency m is contained in f{(f)) or its 
products, demonstrating that the distribution contains products of the frequencies contained 
in /((?!)). Therefore, in the case that the droplet boundary contains only frequency n the 
corresponding state is a superposition of Schur polynomials of dimension kn involving all 
k>0. 

6.3 An example: a disc with a ripple 

Consider a disc with a ripple such that the boundary of the distribution is at 

\w\'^ = 1 + a cos(n<^). (6.25) 

By the arguments above, such a distribution corresponds to a state |$) which is a super- 
position of Schur polynomials of dimension kn, namely 

l*)= E«MA}I^^;W)- (6.26) 

k,{\} 
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In general, one will need to compute all one-point functions of all single and multitrace 
operators in order to deduce the coefficients a^^^xy- However, given that this distribution 
is highly symmetric, with only one defining parameter, let us try to use the vevs of the 
lowest dimension single trace operators to deduce the superposition. Evaluating (3.58) for 
the distribution with boundary at (6.25), the energy and R-charge in this state are 

E-Ec = J=\N'^a^, (6.27) 

whilst from (3.61)-(3.62) the vevs of the neutral operators are 

(c^5-) = ^«^; (05«) = ^«^ (6.28) 



4^/37r2 ' ^ ' 2^/5- 



and from (3.60) the vevs of the maximally charged operators with dimension n are 

Ar2 

{Osn,±n) = ^r-^r^{n - 2)V^rn:eT2mt^_ (g_29) 

All a/ Tl 

First note that the energy and the vevs of neutral operators are reproduced by a superpo- 
sition of the form (6.15) with 

To compute the energy and the neutral vevs at leading order in N we may use any repre- 
sentative orthonormal state \kn), using the result of (6.4). Thus the energy is given by 

2 2 °° f AT ^2fc 

{E - E^U = e-^ J: ^^,{kn) = \nW . (6.31) 

fc=i 

The vev of the dimension four operator can be computed using the result given in (B.6): 



oo 



//o40\ _ (iV«)^S .o40\ .f^ooN 

)^ - e 2n )kn (6.32) 

k=Q 

~ VStt^ ^22W!^''"^-2V57r2" ■ 

Now consider the vev of the maximally charged operator: as we have emphasized this vev is 
sensitive to the specific Schur polynomials contained in the state \kn). Let the single trace 
operator of dimension n be O"^'"" = MnSn, where by construction is orthonormal. Now 

suppose that \kn) is the state created by products of this operator such that 

\kn) = -^{sn)''\n). (6.33) 
V fc! 

The states \kn) are orthonormal to each other in the large N limit. The easiest way to check 
this is to go back to SYM language where \kn) = (^(TrZ")/ViV^) and then use 
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standard large N counting. The leading contribution in {kn\kn) comes from disconnected 
diagrams where each TrZ" is contracted with a corresponding TrZ". 
Then by construction 



{n{k + l)\sn\kn) = VfcTl, (6.34) 

which implies that 

in exact agreement with the holographic result (6.29). By contrast, if \kn) were instead the 
state created by the single trace operator Tr{Z'^"'), then using the appropriate single trace 
operator three point functions one can show that the holographic result for the charged 
operator vev would not be reproduced: the vev would be down by a factor of N. 

Note that when a is infinitesimal the state |$) reduces to a perturbation of the conformal 
vacuum by the operator Sn\^)' 

\^)^{l + ^sn)\n), (6.36) 

in agreement with the identification of infinitesimal perturbations made in [9]. Thus the 
disc with a ripple is consistent with being dual to a state created by coherent superpositions 
of states created by powers of the single trace operators s„ acting on Indeed, the state 
1$) can be written as a coherent state, 

|$) = |5)^e-2^'^4=|fcn), (6.37) 
for S = Naj^^fn with \kvi) the state containing k quanta of s„. 

6.4 Ripplon deformations and the chiral boson description 

The identification of the ripple deformed disc with a state built from a coherent superpo- 
sition of the operators s„ follows naturally from earlier discussions of edge excitations in 
quantum Hall liquids, see for example [39, 40]. Consider a distribution consisting of a single 
droplet, whose boundary can be parametrized as 

= \ + X{^), (6.38) 

where is an arbitrary function with no zero modes. This function X{^^ describes area 
preserving deformations of the disc. Now so far we have used the fermion picture, that is. 
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the Schur polynomials, to describe excitations relative to the conformal vacuum. This is 
a natural basis to describe droplets which are separated from the disc, but it is not the 
natural basis for describing coherent ripplon deformations in the shape of the droplet. 

Such ripplons can best be described by quantizing the chiral boson field X {(/))•, quantizing 
its Fourier mode expansion gives rise to a Hilbert space associated with bosonic creation 
and annihilation operators. As discussed in [40], these operators can in turn be identified 
with elements in the symmetric polynomial (or equivalently, the trace) basis. Thus edge 
waves and deformations in the shape of the droplet are most naturally described within the 
trace (chiral boson) description rather than the Schur polynomial (fermion) description. 
This extends the identification made for the single frequency ripple discussed above to more 
general ripples. 

It is interesting to note that the algebra of area preserving diffeomorphisms of the droplet 
is actually Woo- It emerges in both the fermionic and bosonic formulations as the algebra 
of unitary transformations of physical states. It would be interesting to understand the 
meaning and implications of W^o for holography. 

7 Discussion 

In this paper wc have discussed holography for bubbling solutions. Sohitions that are asymp- 
totically AdSp X S** contain an infinite amount of holographic data that can be extracted 
by algebraic manipulations, namely the holographic 1-point functions that characterize the 
vacuum of the dual QFT. This is the simplest information one can extract from a given 
supergravity solution. Conversely, knowledge of the 1-point functions is in principle suffi- 
cient in order to reconstruct bulk solutions from QFT data. Two-point and higher-point 
functions can also in principle be extracted, but this requires solving at least the linearized 
equations around the solution (for 2-point functions) whilst for n-point functions one needs 
to solve the (n— l)-th order equations. Explicitly solving such equations is an intractable 
problem, except for very symmetric solutions. Moreover, in the case of interest vevs are 
protected, given the non-renormalization of 3-point functions of chiral primaries at the con- 
formal vacuum, whilst there are no such non-renormalization theorems protecting generic 
two point functions, and corresponding four point functions in the conformal vacuum. 

In the first part of this paper we reviewed the holographic 1-point functions derived in 
[4] for the stress energy tensor, the R-currcnts and all chiral primaries up to dimension four 
for asymptotically AdS^ x solutions of IIB supergravity that involve the metric and the 
5-form. These results hold generally when the solution is dual to a state (rather than describ- 
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ing a deformation), i.e. they do not dependent on the amount of supersymmetry preserved 
by the solution or its bosonic isometries (except that the solution should be asymptotic 
to AdSs X S^), and thus these 1-point functions can be used to extract holographic data 
from any such solution. The 1-point functions are given in terms of the asymptotic coef- 
ficients of the ten dimensional solution and are presented in (2.25), (2.35) (2.36) (which 
for the reader's convenience were boxed). Note that the expressions are non-linear in the 
asymptotic coefficients. 

The next step was the evaluation of these holographic formulas on the LLM solutions. 
The holographic 1-point functions are by construction diffeomorphism covariant, so the 
asymptotic coefficients can be extracted in any coordinate system. A clever choice of 
coordinates however can reduce the required labor significantly. Recall that the LLM 
solution is determined by a harmonic function $ in six dimensions with sources on a 2- 
plane. The asymptotic expansion around AdSr^ x can be efficiently performed by writing 
$ = $° + A$, where $° is the harmonic function that leads to AdS^ x S'\ and then ex- 
panding in A#. One still needs to convert these expansions into radial expansions. This is 
done by first using fiat coordinates on so that the asymptotic expansion of A$ takes a 
standard form and then transforming to the coordinates most natural for the AdS^ x so- 
lution. This procedure minimizes the number of non-linear terms entering the computation 
of the 1-point functions. 

We obtained explicit expressions for the vevs in terms of integrals over the 2-plane 
defining the solution; these are given in (3.56)-(3.58)-(3.60)-(3.61)-(3.62) (which again are 
boxed). Note that in (3.60) we give the vevs of all maximally charged operators, i.e. op- 
erators with dimension equal to R-charge, despite the fact that the general analysis in the 
previous section was done for operators up to dimension four. The ability to produce such a 
result is due to special properties of the LLM solution combined with the previous Coulomb 
branch results of [2]. Given the large amount of supersymmetry preserved by the LLM 
solutions, one would expect that these vevs should not rcnormalizc and thus that they must 
be reproduced by a weak coupling computation. Put differently, these vevs provide checks 
for the correct identification of the dual theory. 

The vevs satisfy a number of non-trivial consistency checks. Firstly, the vcv of the 
energy is proportional to the vev of the R-charge (up to the the Casimir energy of SYM 
on 5^) as is required by supersymmetry. Secondly, all vevs, except for the energy which 
should become equal to the Casimir energy, should vanish for the theory at the conformal 
vacuum and this is indeed the case for the vevs we derive. Thirdly, the LLM solutions in the 
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decompactification limit of go over to SO{4) symmetric distributions of D3 branes, with 
the sources on the 2-plane now describing the distribution of D3 branes. These solutions 
are dual to N = 4 SYM on the Coulomb branch. The LLM vevs indeed correctly reduce to 
the Coulomb branch vevs given in [2] in this limit. 

Before proceeding we should make a comment about the mass of these spacetimes^. 
Had the solutions been asymptotically fiat, one would have obtained their mass from the 
gtt component of the metric. In our case however the solution is asymptotically AdS^ x 
so such a prescription is in general not valid. There are two issues here. One is that the 
solution involves in a non-trivial fashion a compact part of the geometry and the other is 
that the non-compact part is asymptotically AdS. For asymptotically AdS spacetimes the 
issue of mass has been revisited and thoroughly analyzed in recent years [23, 24, 25, 44, 45] 
resulting in holographic formulas which relate the mass to the asymptotics of the metric 
and other matter fields. 

Taking into account the compact part is also non-trivial since none of the existing 
consistent truncation formulas from ten to five dimensions is directly applicable. One can 
however reduce the solution to five dimensions without truncating, keeping all fields relevant 
for the computation of the mass. This is essentially the method of KK holography [4] and 
results in the rigorous formula for the holographic stress energy tensor given in (2.35). One 
can then obtain the mass from the Tu component, as usual. 

The field theory dual to the LLM solutions is expected to be iV = 4 SYM on R x 
in a half supersymmetric state. A general way of analyzing this theory would be to carry 
out path integral quantization. The requisite supersymmetry is preserved by quantizing 
around 1/2 supersymmetric solutions of AT = 4 SYM on Rx S^. Examples of such solutions 
were discussed in [42] and more recently in [43]. These solutions are time-dependent and 
in correspondence with the Coulomb branch of AT = 4 SYM on R^^'^\ In particular, the 
curvature coupling implies the scalar Z satisfies an equation of the form Z ^ iZ. This 
implies that the R-symmctry current ~ TrZ Z and the operator O'^'^ ~ TvZZ 
are proportional to the each other when evaluated on these solutions (and similarly for 
related higher dimension operators). This provides an additional consistency check for the 
holographic vevs, which the vevs in (3.56)-(3.61) indeed satisfy. 

Due to the extended supersymmetry one might expect that the exact values of the 
vevs of 1/2 supersymmetric gauge invariant operators could be computed by a semiclassical 
computation. This would provide a rigorous computation of the vevs from first princi- 
*An earlier discussion about the mass of the LLM solutions can be found in [41]. 
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pies. Actually carrying out this computation is not so easy in practice, though, because 
of subtleties associated with correctly treating the integration measure. The issue is the 
following: for the computations of interest one will want to integrate out most of the SYM 
fields, including the other four scalars X* and off-diagonal degrees of freedom of the complex 
matrix Z, leaving only an integral over the eigenvalues of this matrix. This in turn involves 
correctly parametrizing the path integral measure as given in 



and then integrating out appropriately. Now for a general computation one does not expect 
to be able to integrate out exactly all these degrees of freedom. Integrating first over the 
S"^ would lead to a complicated interacting multi-matrix model which will not in general 
be solvable. However, the holographic computations for the vcvs along with the fact that 
we can reproduce them by the holomorphic matrix model, imply that at least for these 
computations one can explicitly integrate out the other degrees of freedom. Demonstrating 
this by a first principles computation would be Tiscful as it would explain the regime of 
validity and the limitations of the free fermion description. Moreover, one may in this way 
show how certain computations can be carried out in multi-matrix models, even when they 
are not exactly solvable. 

In the absence of a rigorous derivation of the free fermion description from first principles, 
we proceeded by using it as a working assumption. On symmetry grounds the state that any 
given bubbling solution is dual to is a superposition of states obtained from the conformal 
vacuum by the action of a 1/2 BPS operator. The question is then whether one can 
uniquely determine the precise superposition from the data encoded in the solution. Using 
the identification of the coloring of the 2-plane with the phase space distribution of the free 
fermions we show that this information alone does not completely determine the state in 
the large N limit. It does determine it enough however so that the vevs of all single trace 
1/2 BPS operators in that state are uniquely determined. This is precisely the information 
encoded holographically in the asymptotics of the solution. The missing information is 
related to vcvs of multi-trace operators. 

A general single trace 1/2 BPS operator depends on fields other than the complex Z 
field. This implies that these operators cannot not in general be implemented with free 
fermions. Nevertheless, we showed that for the purpose of the computation of the vevs 
and to leading order in such an implementation is possible and all such operators are 
expressed in terms of bilinears of fermion creation and annihilation operators. Using these 
expressions we show that all vevs computed holographically agree exactly with the field 




(7.1) 
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theory computation in the large N hmit for any distribution. 

To illustrate our discussion we analyzed a number of examples. In accordance with our 
general discussion, we showed that all vevs associated with any symmetric distributions are 
degenerate to leading order N. For non-symmetric distributions, the vevs of charged oper- 
ators (which by symmetry considerations are zero in symmetric distributions) can (partly) 
distinguish between different states. However, an infinite superposition of states of definite 
R-charges is required to obtain a regular geometry. We also analyzed in detail the case of 
the distribution being a ripple on a disc. This case has been analyzed previously for an 
infinitesimal ripple in [9]. We showed here that a finite ripple corresponds to a coherent 
state of single trace operators. 

We should also comment on the striking parallels between this system and the 2-charge 
D1-D5 fuzzball solutions. Both systems can be characterized by a set of curves: in the LLM 
case these are curves in describing the droplet boundaries, whilst in the D1-D5 case these 
are curves in an auxiliary space describing the supertube shape and its internal excitations. 
The holographic analysis for this system has recently been done in [46]. In both cases, only 
when the curves are circular and preserve rotational symmetry do the geometries correspond 
to vacua built from a single operator (in the R-charge basis). Regular geometries in which 
the rotational symmetry is broken correspond to infinite superpositions of states in the 
R-charge basis, with the coefHcients of the superpositions related to the Fourier expansions 
of the curves. Thus the natural bases in the dual field theory, which are labeled by their R 
charges, are not the natural bases for regular geometries. 

It would be interesting to use the holographic anatomy techniques discussed in this 
paper to analyze 1/4 and 1/8 BPS bubbling solutions [34]. One would expect that these 
include both geometries dual to states and those dual to deformations. A holographic anal- 
ysis should determine how the boundary conditions which ensure regularity in the interior 
of these geometries are related to the vevs/deformations in the dual theory. More generally 
one may hope that combining supcrsymmctric classification techniques with holographic 
anatomy might lead to more efficient holographic engineering of geometries dual to super- 
symmetric field theory states and deformations. 
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A Properties of spherical harmonics 

The defining equations for the spherical harmonics are 

UyY^^ = A^^Y^\ A^^ = -k{k + 4), fe = 0,1,2,... (A.l) 

= ^hyh^ A'^ = -{k^ + 4k-l), k = l,2,... 
Dy^S = ^'''Y^abV A'^^=-(fc' + 4fc-2), A; = 2,3,... 

^ ^'"^fi' A'- = -(A:^ + 4fc-2), A; = 1,2,... 
OTi^ = OT^^^i,'; = = 0. 

The overall normalization is chosen so that the harmonics are normalized as 

The triple overlap between spherical harmonics is defined as 

J y(fei,/i)y(fe2,/2)y(fe3,/3) = n^ak,l,MIM, (A.3) 

where {k, I) is the degree of the scalar harmonic and / labels the SO{6) quantum numbers. 
Recall that the scalar harmonics can be represented as 

yikuh) ^ (jh x'^...x^^ (A.4) 

*l"''tfc ^ ' 

where x'" are Cartesian coordinates on and Cf^,,,^^ is a totally symmetric traceless rank 
k tensor of 50(6). The normalization in (A.2) corresponds to delta function normalization 
for the C^'s, i.e. 

Note that 

(^S + 2)!22(^~2) ai!a2!a3! 
where T, = ki + k2 + kz-, 0:1=5(^2 + kz — ki) etc. Useful identities for the scalar harmonics 
include 

D''D(,D^)Y' = 4{1 + ^)DaY'; (A.7) 
nyD^,D,)Y' = {W + A')D^,D,)Y'; 
DyDaY^ = {A^ + 4)DaY^. 
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Vector harmonics are normalized so that 

j Y^^Y^^" = TT^z{k)5^^^\ (A.8) 

where z{k) is as given in (A. 2). We introduce the fohowing coordinates on 

ds^ = d9^ + cos^ 9dnl + sin^ edcf)'^. (A.9) 

The differential equation (A.l) for the scalar harmonics is separable. Imposing 50(4) 

symmetry implies that the spherical harmonics depend only on 6 and 4>. The general 
solution can then be expressed in terms of a hypergeometric functions, 

y('='-)(e,<^) = C(„,^)y^(e)e^-'^ (A.IO) 

where C(^n,m) ^ normalization constant and the function y^{0) is given by 



{x) = a;l"^liF2(-^(fc - \m\), 2 + ^(A; + |m|), 1 + \m\-x^) (A.ll) 



with X = smO (there are also a second solution with leading behavior a;^'™' but this solution 
docs not reduces to a finite polynomial for any choice of the quantum numbers). The 
hypergeometric function reduces to a finite polynomial when cither the first or second 
argument is zero or a negative integer. This leads to the following cases 

(A; = 2/, m = 2n), {k = 2l + l, m = 2n+l) n G Z G Z+ (A.12) 

with 

yii(x) = x2WiF2(-/ + |n|,2 + / + |n|,l + 2|n|;x2) (A.13) 

yln+ii'') = a;|2"+iliF2(-/ + |n|,3 + / + |n|,2 + 2|n|;x2^ 



The harmonics that are also SO (2) symmetric are given by 



The lowest harmonics are therefore 



y('''^ = -l=(3sin2 0-l), (A.15) 

y(4,o) ^ ^(10sin^^-8sin2e + 1), 
4v5 
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We will also need the following normalized charged scalar harmonics 

y(fe.±*^) = J_ sin*^ ^e^^'^'^; (A.16) 

y(3,±i) ^ ^sin0(2sin2^-l)e±^'^; (A.17) 

y{4,±2) ^ -^sm^e(5sm^9-3)e^^'^: (A.18) 
2710 



Note that the triple overlap between charged and neutral harmonics is given by 

,fJ{k,-k)(J{k,k)(Ji2p,0)^ ^ ^ (A.19) 

where C^^''^^ denotes the symmetric tensor corresponding to the degree p, SO{2) charge q 
spherical harmonic. The relevant vector harmonics are those with only components along 



= ^sin^^d^; (A.20) 
v2 

/o 

= ^sm^e{2sm^e-l)d(l>. (A.21) 
In extracting the perturbations the following are useful: 

D^eDo^Y^'^^ = -i=e±2^<^(^-48sin4^ + ^sin2 0-6^; (A.22) 

D^^D^^Y^'^-" = ^e±2'^sin2^('l2sin^e-^sin2e + 6 



and 



1 


2^ 


/lO 




1 


2^ 


/lO 




1 



5 

Dr^D^.Y^'^^ = ^e±2^'^cos2dl2sin4^-?sin2| 



2V10 V 5 

DeD^Y^^"^^ = ±i-^e^2#(15sin3ecos0-3sin^cose), 
v 10 



D^eDe)Y^'^' = e^^^^(^-Um'e + iy, (A.23) 
DuD^-,Y^'^^ = e^^^"^ sin^ ^ ( - sin2 ^ - 1 



De-D^y^'^^ = ize^^^'^sin^cos^, 
along with the corresponding expression for the neutral harmonics: 

D^eDe)Y'' = ^(-96sin^e + ^sin2^-|); (A.24) 



An If 4 24 o 48 

D,,D,,Y^^ = ^(24sin^^-ll^sin^^ + | 
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B Scalar chiral primaries 



The single trace scalar chiral primary operators of dimension k are defined as 

where the properties of the degree k symmetric traceless tensors Cj^,,,-^ are given in appendix 
A. The operators are normalized such that 

where the scalar fields X*" are normalized such that 



{X^{x)X^{y)) = (B.3) 

y\ 

where (a, 6) are color indices. The appropriate normalization of the dimension k chiral 
primaries to match with supergravity is 

^fi = -^{k-l)ik-2f, (B.4) 

for A; 7^ 2 with Af^ = A^^/ttI 

The planar three point function for such scalar chiral primaries is given by 

Here 2a3 = ki + k2 — k-^ and (ai, 02) arc defined analogously. The triple overlap of the sym- 
metric traceless tensors is denoted {C^^C^^C^^}; recall that these tensors are orthonormal 
(A.5). 

Now let us consider the specific case of three point functions between one neutral 
(S'0(4) X 50(2) singlet) operator O^fe.o ^j^j^ dimension 2k and two conjugate S0{2) charged 
operators O"'". The corresponding spherical harmonics are given in (A. 14) and (A. 16) re- 
spectively, with the triple overlap being given in (A. 19). The three point function implies 
that the vev of the neutral operator in the (unit normalized) state created by O"'" is 



Therefore the vevs of the neutral operators in these states are given by 



/n^fi\ _ . /n^k,o^ _ n {k-1) / 2/g(2fc- 1) 
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We will make use of several other three point functions, involving two maximally charged 
operators: 

(03,-3^3,102,2^ ^ :^3V2(C(3'-3)C(3.1)C(2.2)) = V3^; (B.8) 

Ma . rz I „C4 -A\ „C4 . A/4 



(04,-4(^4,2(^2,2^ ^ :^4^((^(4,-4)^(4,2)^(2,2)^ ^ 4. 



where C^'^) denotes the symmetric tensor corresponding to the degree p, 50(2) charge q 
spherical harmonic. 



C Large N behavior of three point functions 

To compute vevs of single trace chiral primary operators in generic half BPS states we use 
the corresponding three point functions. To determine the dominant effects in the large N 
limit we thus need to know the N dependence of correlators of the form 

C,:C,jJ{Tr:{aiZ)rO^Tr:{aiZ)) (C.l) 

where O"^ is a single trace chiral primary, denotes a conjugacy class of Sn and the 
normalization factors C^i^ are such that the operators are orthonormal in the large N limit: 

C,,C,jJ{T,{aiZ)rT.{aiZ)) = SnmS'' + 0{l/N). (C.2) 

Using the propagators given in (B.3) one finds that C^j ~ l/N"'; note that throughout this 
section we will suppress factors of order one. It is convenient to introduce the notation 

i i i 

for an operator of dimension n involving m traces with a permutation labeled by a[m\. 
As discussed in [36] there are two distinct cases of correlators to consider, the extremal 
correlators in which the dimension of the conjugate operator is equal to the sum of the 
dimensions of the other operators and non-extremal correlators. 

Let us consider first non-extremal correlators, focusing on the case where O"^ is an 
SO{2) neutral operator, namely it is an operator O'^^'^ of dimension 2p such that 

O^P''' = ^Tr{ZPZP + ■■■), (C.4) 

where the ellipses denote cyclic permutations. Now charge conservation implies that the 
correlator 

(^0a[mi]-)t02p,O(0a[m2])^ (C.5) 
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is only non-zero when ni = n2- The N dependence varies according to the specific choices 
of ((T[mi], cr[m2]). As discussed in [36], for a generic choice the correlator will have the same 
N scaling as the related (mi + m2 + l)-point correlator of single trace operators, namely 
as l/_/V'"i+'"2-i. ^]jug foj. single trace operators the scaling is 1/N. Recall that an n-point 
correlator of single trace operators behaves as 

^QklQk2 ...Qkn^ ^ __L^^ „>2. (C.6) 

However, for specific choices of (cr[mi], cj[m2]) the N scaling can be enhanced, because there 
are disconnected components to the diagrams. In particular, large counting gives 

^(0aM)t02p,O(0aH)^ ~ ^, (C.7) 

for any m, whilst for (T[mi] / cr[m2] one always finds a subleading N dependence, with the 
3-point function being at most of order 

(This result for m\ = l was given in [36] .) Thus vevs of neutral operators are thus dominated 
by diagonal three point functions of the type (C.7). For the vevs of non- maximally charged 
operators, the relevant correlators are also non-extremal; the leading terms scale as XjN and 
arise from single trace correlators and specific multi-trace correlators. We will not however 
need detailed results for the latter. 

Now let us turn to the extremal correlators in which O"^ is a maximally charged single 
trace operator. Again the correlator involving single trace operators behaves as l/AT, but 
in this case correlators involving multi trace operators can dominate, since they can grow 
as 1. In particular, 

((cS^'')^C''''*^(Cn^"'^'')) ~ 1, C'^il^'^ = C»*^''=(C»^H). (C.9) 

Note that analogous results are obtained in the Schur polynomial basis; see [6] for related 
discussions. 

D Killing spinors for LLM solutions 

We discuss in this appendix the computation of the Killing spinors of the LLM solutions. 
This computation was carried out in appendix A of [5] but only half of Killing spinors 
were correctly identified, even though the projection operators were given correctly, and 
furthermore the spacetime dependence is not given correctly. These corrections do not 
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affect tlie final answer for the supergravity solution (although some intermediate steps in 
the derivation are affected). They may have a real effect however in similar computations for 
less supersymmetric solutions. Furthermore the correct Killing spinors may be needed for 
other purposes, for example for analyzing supersymmetric probe branes in this background. 
We use the notation of [5] and choose the same basis of gamma matrices 

r^u = 7m ® 1 <^ 1 <^ 1' r„ = 75 (g) cr„ (g) 1 (g) ai, = 75 (g) 1 (g) a„ (g ^2, (D.l) 

where aa,d'a, (^a arc the Pauli matrices. 

The ten dimensional spinor is decomposed as 

where Xa-, Xa geometric Killing spinors of S^, i.e. they obey 

^cXa = a-^lcXa, a = ±1, (D-3) 

and a similar equation for %„, where Vc is the standard connection on a unit 3-sphere. We 
normalize these spinors as XaXa = XaXa = 1- The fact that the spinors are correlated as in 
(D.2) follows from the analysis in [5]. 

The Killing spinor equation then reduces to [5] 

(zae-5(^+G)^5^^ + lJ^'^^,{H + G))e + 2Me = 0, (D.4) 

(zae~2(^-G)75CT2 + \l^d^{H - G))e - 2Me = 0, (D.5) 

V^e + M7^€ = (D.6) 

where 

M = -^e-i(^+«)7'^''F^,75ai (D.7) 

Processing these equations one finds that the spinor should satisfy the following equations 
[5] 

p-ea = Rtea = 0, (D.8) 
where we introduce the commuting projection operators 

P± = i (1 ± (ie-^75 + aVl + e-2Gr3<Ti)) , = i(l ± mFiFs) (D.9) 

satisfying 
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Each of this projection cuts the number of spinors by 1/2, so we have a total of 8 Kilhng 
spinors for a = +1 and 8 Kilhng spinors for o = —1. The most general solution of (D.8) is 

e„ = R-P^ia (D.ll) 

where ia are (at this point) unconstrained spinors. 
In [5] the following solution of (D.8) was given, 

e = e^'^^'^'^'ei, rVei = aei, sinh2,5 = ae-^ (D.12) 

These are in fact only half of the Killing spinors in (D.ll). To see this introduce a new 
projector, 

5^ = i(l±arV), {Stf = St StS- = 0, (D.13) 
and decompose as 

e„ = e+ + 6-, S±6t = et, Ste^ = (D.14) 
A short computation yields, 

P+e+ = cosh(5e^'5T'r''^'e+ (D.15) 

which is the spinor in (D.12). Upon multiplication by R~ one obtains half of the Killing 
spinors in (D.ll), namely we miss the ones based on e~. 

To specify the Killing spinor we need to specify To this end we consider the fermion 
bilinear /2 = iea2e. It was shown in [5] that /2 equals 

Inserting the spinors in (D.ll) and defining 

~et = c^et (D.17) 



we find that (D.16) implies 



and 



where 



c± = - (D.18) 

VVl + e-2G ± 1 



ieaa2R~ea = 2 (D.19) 
e„ = e+ + n5e-, (D.20) 
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The Killing spinor becomes 

Ca = — ^e4(^+'^)i?~ ((cosh + ia75 sinh(5)e^ + (— sinh(5 + ia75 cosh 5)e~) 
v2 

= ^ei(^+«)i?- {e^'^'^'^'ei + iaj^e'^'^'^'^'e-) (D.21) 

Prom these spinors one can construct appropriate fermion bilinears and determine the su- 
pergravity solution as in [5]; the supergravity solution is exactly as given in [5]. Note that 
the functions {G, H) are such that 

e^ = y\ z = itanhG, (D.22) 

where z is the defining function of the supergravity solution. 

There is however a further issue in constructing the actual Killing spinors: the spinors by 
construction satisfy (D.4) and (D.5) since it is these equations which were processed. One 
still needs to check explicitly that all components of (D.6) are satisfied. Now the spinors 
as given in [5] do not depend at all on the time coordinate t. This is however inconsistent 
with the t component of (D.6); one can show that 

Vte„ + M7te„ ^ (D.23) 

for constant (e+ , e~ ) . Another way to see that the Killing spinor solution is not quite correct 
is by considering the limiting case of AdS^ x S^. The known explicit expressions for the 
Killing spinors of AdS^ x do depend explicitly on the time coordinate; this remains true 
for the specific combinations of spinors which form the set of sixteen discussed above. 

The resolution of this issue is straightforward: {c-t^'^a) not constant, but must 
contain suitable t (and indeed also (p) dependent phase factors so that (D.6) is satisfied. 
These phase factors drop out of (D.19) and all other fermion bilinears used to construct the 
supergravity solution. 
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